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1- Tempo entre Falhas 


Os dados de falha fornecidos para o equipamento hipotético encontravam-se desordenados, de forma que, 
para uma avaliação estatística foi necessário tratá-los, ordenando-os em ordem crescente, e expressando-os em 
horas decimais. Para tal, foi calculado o tempo decorrido entre uma falha e outra, a partir do início do 


funcionamento do maquinário. Os valores obtidos estão na tabela 1, abaixo. Para efeito de cálculo, foi utilizado 


o Excel. 


A manutenção preditiva é uma manutenção planejada que na qual através de coleta de dados e medições, 
feitas em um intervalo regular de tempo, pode-se conhecer as condições de cada um dos componentes das 
máquinas e equipamentos, sendo assim podese ser feito o acompanhamento do desgaste de peças vitais de um 
equipamento, assim como sua vida útil. Esta manutenção é a de menor custo. Os dados abaixo representam dados 


de Falhas e Reparos de um equipamento hipotético, que podemos analisar para subsidiar na elaboração de um 


plano de manutenção preventiva do referido equipamento. 


Dados do Aluno Arthur Sena e Silva - 20190070675 


Tabela 1 - Dados do aluno Arthur Sena 


Tempos de Reparos Hora da Falha 
[Hora:Minuto] Datas das Falhas | [Hora:Minuto] 
14:44 11/01/2020 07:40 
09:20 20/01/2020 23:25 
10:45 28/01/2020 00:27 
13:50 02/02/2020 1511 
10:24 11/02/2020 20:15 
18:42 17/02/2020 18:32 
04:05 26/02/2020 03:43 
04:32 05/03/2020 03:41 
07:34 15/03/2020 04:33 
10:02 25/03/2020 10:09 
14:50 01/04/2020 [ON ZON 
05:40 08/04/2020 04:57 
10:46 13/04/2020 21:43 
10:08 23/04/2020 08:24 
23:56 28/04/2020 23:58 
03:50 05/05/2020 2121 
10:36 14/05/2020 00:59 
07:25 22/05/2020 18:47 
12:13 31/05/2020 06:17 
04:43 07/06/2020 15:49 
59:35 15/06/2020 00:14 


05:32 27/06/2020 05:37 
09:11 07/07/2020 19:28 
LEA 14/07/2020 20:52 
06:49 23/07/2020 11:14 
05:25 02/08/2020 18:50 
1353 11/08/2020 12:54 
09:22 21/08/2020 02:41 
05:57 26/08/2020 14:20 
03:54 04/09/2020 07:24 
LH SA 14/09/2020 02:18 
07:53 22/09/2020 07:10 
04:36 30/09/2020 15:48 
02:15 07/10/2020 08:44 
08:21 16/10/2020 09:44 


1.1 AJUSTE DOS DADOS E CÁLCULO DOS TEMPOS PARA FALHA 


O presente trabalho tem como objetivo analisar e dar subsídios na elaboração de um plano de manutenção 
preventiva de um equipamento hipotético. Inicialmente os dados de falha foram reorganizados, 
posteriormente foram calculados os tempos para falha e tempos de reparo. Os resultados obtidos são 
mostrados na tabela 1. Para encontrar o primeiro tempo para falha considerou-se que o equipamento entrou 


em funcionamento no dia 01/01/20120, às 8:00 horas. 
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Tabela 2 - dados de média, moda/ mediana e variância 


683.95 11:30:38 
718.65 08:34:00 
x y 
14902.61 05:09:57 
123.91 11:08:05 


2. Métodos Empíricos. 


Os métodos empíricos foram utilizados para uma estimativa inicial dos parâmetros das distribuições tanto 
de tempos de falhas quanto nos tempos de reparos. Como primeira hipótese foi suposto que os 35 dados não 


sofreram censura e, dessa forma os rank”s médio e mediano puderam ser utilizados. 
2.1. Falha pelo método Rank Médio de Herd-Johnson. 


Usando o método rank Médio de Herd-Johnson, técnica estatística utilizada para analisar dados de 
classificação. Ele foi desenvolvido por John Herd e Alan Johnson é particularmente útil em situações onde os 
dados estão sujeitos a erros ou imprecisões. Além de ser aplicado quando os dados não estão precisamente 
classificados em ordem, mas sim em grupos ou faixas. Em vez de calcular a média dos valores numéricos em si, 
como no método da média aritmética, o método Rank Médio calcula a média dos postos ou classificações dos 


valores 


O rank médio é comumente utilizado para estimar o desenvolvimento da função cumulativa de 


falha F(t) que é função apenas da posição do dado depois de ordenados de forma crescente, expressa por: 


O rank médio é comumente utilizado para estimar o desenvolvimento da função cumulativa de falha F(t) que é 
função apenas da posição do dado depois de ordenados de forma crescente, expressa por: 
i 


Pt)=1-R(t)= n+1 


Sendo n o número de falhas no intervalo observado e i a posição da falha nos dados ordenados. A 
confiabilidade, portanto, é dada por: 


i on+1-1 


Rt)=1-— = 
(ti) n+1 n+1 


Para a caracterização da distribuição de probabilidade foi estimada a função densidade de probabilidade 


dada pela equação: 


1 


ft) = AG +1) 


Com isso obtém-se a taxa de falha, expressa por: 


Sê) 1 
AG = RPC) Atin+1-i) 


O tempo médio entre falhas (TMEF) foi estimado através da equação 


t; 
TMEF = a = 717,32 horas 


O tempo médio para reparo (TMPR) foi calculado a partir da estimativa da taxa de reparo: 


K 


L=57— 0,171 
Onde: 
K: Número de reparos; 
tr: tempos de reparo; 
Daí: 
TMPR = E =-——— = 10,17 horas 
u 0,1051 , 


Disponibilidade = TMEF/ (TMEF + TMPR) 
= 683,95/(683,95 + 11,5) = 0,1017 horas 


Os valores obtidos pelas estimativas são mostrados na tabela 2 e as figuras de 1 a 4. 


Tabela 3 - Utilização do Rank Médio para obtenção da função cumulativa, confiabilidade, taxa de falha e densidade de falha, através dos tempos 
entre falha. 


11 129,2090278 129,2090278 47.00208 0.305556 0.694444 0.000851 0.000591 
12 138,5395833 1385395833 0.173611 0.333 0.666667 0.24 0.16 
13 86,13819444 86,13819444 0.151389 0.361111 0.638889 0.287196 0.183486 
14 206,2166667 206,2166667 0.354167 0.388889 0.611111 0.128342 0.078431 
15 73,06527778 73,06578 23.41389 0.416667 0.583333 0.002034 0.001186 
16 14105625 14105625 0.152083 0.444 0.555556 0.328767 0.182648 
17 181,4409722 181,4409722 22.46597 0.472 0.527778 0.002343 0.001236 
18 161,3659722 161,3659722 0.660417 0.5 0.5 0.084122 0.042061 
19 180,0451389 180,0451389 0.652083 0.5278 0.472222 0.090209 0.042599 
20 139,9423611 139,9423611 0.240278 0.5556 0.444444 0.260116 0.115607 
21 156,4263889 156,4263889 0.529861 0.5833 0.416667 0.125819 0.052425 
22 258,7006944 258,7006944 26.16597 0.6111 0.388889 0.00273 0.001062 
23 207,6277778 207,6278 21.60625 0.6389 0.361111 0.00356 0.001286 
24 1318527778 1318527778 0.308333 0.666667 0.333333 0.27027 0.09009 
25 185,6513889 185,6513889 0.353472 0.694444 0.305556 0.257189 0.078585 
26 206,7451389 791.2951 22.45417 0.722222 0.277778 0.004454 0.001237 
27 90,64722222 792.0743 0.779167 0.75 0.25 0.142602 0.035651 
28 188,4083333 814.9854 22.91111 0.7778 0.2222 0.005456 0.001212 
29 203,2291667 816.1035 1.118056 0.80556 0.194444 0.127773 0.024845 
30 160,4152778 816.3097 0.20625 0.833 0.166667 0.808081 0.13468 
31 164,0583333 839.9076 23.59792 0.8611 0.138889 0.008475 0.001177 
32 1421138889 839.9132 0.005556 0.888889 0.111111 45 5 

33 206,7451389 840.1826 0.269444 0.916667 0.083333 1.237113 0.103093 
34 90,64722222 840.4889 0.30625 0.944444 0.055556 1.632653 0.090703 
35 840.49 0 0.972222 0.027778 


Abaixo, é mostrado como foi o comportamento da função cumulativa, pelo Rank Médio dos dados de 


falha. 
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Figura 3 - Comportamento da Função Cumulativa, pelo Rank Médio dos dados de falha. 
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Os valores obtidos pelas estimativas pelo método de falhas são mostrados na tabela 2 e as figuras de 1 a 
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Figura 5 - FIGURA 1: 
Comportamento da Função 
Cumulativa, pelo Rank Médio 
dos dados de falha. 
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Figura 7 - Comportamento da Confiabilidade, pelo Rank Médio dos dados de falha. 
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Figura 6 - Comportamento da Taxa de Falha, pelo Rank Médio dos dados de falha 
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Tabela 3 - Valores obtidos através dos Métodos Empíricos e do Rank Mediano. 


2.2 MÉTODO DE ESTIMAÇÃO USANDO A TÉCNICA DE MÁXIMA VEROSSIMILHANÇA 


Esta seção trará os mesmos resultados da anterior, entretanto os valores serão obtidos com auxílio do 
software ProconfQ pela técnica de máxima verossimilhança. Dois testes de aderência foram utilizados, o qui- 


quadrado e o de Kolmogarov-Smirnov, sendo estes os mais difundidos para aplicação. 


O método da Verossimilhança consiste em estimar parâmetros que permitam a associação entre a 
distribuição dos dados e as distribuições estatísticas padronizadas. Nesse sentido, um vetor de parâmetros L é 


calculado, de tal forma que a aproximação a seguir seja possível: 
L(61,02,03,04...0n)=[|f(ti/01,02,603,04...0n) 


Para cada distribuição em particular, existem equações para a estimação dos parâmetros 01, 02, 63, 04 
... On. Como a resolução analítica é laboriosa, foi utilizado o software ProConf. O programa é capaz de 
fornecer instantaneamente os parâmetros de verossimilhança entre a distribuição da amostra e as distribuições 
comuns, especificamente, Exponencial, Weibull, Gamma, Normal e Lognormal, bem como realizar testes de 
aderência, como será detalhado no tópico seguinte. A primeira etapa foi aquisição dos dados ordenados no Excel, 


conforme a figura 5 


Ao inserir os dados de falha no programa, é possível identificar quais métodos de distribuição se adequam 
ao caso em análise. A tabela a seguir mostra o nível de significância de cada distribuição e, também, se é possível 


aplicá-la ou não. 


Modelos | Ajuste/E statísticas | Funções de Confiabilidade | Gráficos Testes de Aderência | 


. | Tempo até a Falha | 4) 
Detectar o tamanho da amostra e 
E Desa classificar dados em ordem crescente - Processar | | 


Tamanho da amostra Amostra completa 
Número de falhas Censurada por nr. unidades 
Censurada por tempo 
Censura aleatória 


Teste Analítico do 72: 


Limite Inferior | Limite Superior re po Eópadi do 
411145 

561075 
561075 


711005 


Parâmetro de localização 711005 
(* Fomecido pelo usuário Par, localização 78597. 


€ Fomecido pelo computador 0 Z2= 86,48 com 4 graus de liberdade Nível de Significância = 0,0001 


; Teste Analítico de Kolmogorov-Smirnov (K-5): 
— Simular. | DN=0,4681 Nível de Significância = 0,0001 


tra : A hipótese de que a população segue o modelo 
Tecle 'c' para indicar tempo completo ou tempo censurado E RES 
Tecle 'p' para reproduzir o valor da célula anterior Exponencial é rejeitada. 


[Ei] &) Bj Re 


Figura 9 - método wiebull 


Modelos | Ajuste/Estatísticas | Funções de Confiabilidade | Gráficos Testes de Aderência | jFigura 10 - método gama 


Teste Analítico do X2: 


Frequência Frequência Valor do 


Limite Inferior | Limite Superior bssrvada Estao srietiedo 


411145 


Z2= 4,90 com 3 graus de liberdade Nível de Significância = 0,1793 


Teste Analítico de Kolmogorov-Smirnov (K-S): 
DN = 0,1725 Nível de Significância = 0,0099 


A hipótese de que a população segue o modelo Gamma é 
rejeitada. 


[E] &) aj Re, 


Modelos  Ajuste/E statísticas | Funções de Confiabilidade | Gráficos | Testes de Aderência | 


HO =52216.713 
t50 = 67591,634 

MTTF = 68432.632 
Std Dev. = 13153,193 


Distribuição: Gamma 


Parâmetro de Localização = O 


Y= 27.0685 
Var(y)= 58.2576 


95% do Intervalo de Confiança: 
para y= 16.615 até 43.2061 


6= 2528.128 
Var(6)= 353079,746 


Teste Analítico do 72: 
Limite Inferior | Limite Superior mst iso En 
411145 
561075 
561075 
711005 


%2= 5.69 com 3 graus de liberdade Nível de Significância = 0,1279 
Teste Analítico de Kolmogorov-Smirnov (K-S): 
DN = 0,1842 Nível de Significância = 0,0041 


A hipótese de que a população segue o modelo Lognormal 
é rejeitada. 


Analisando o modelo Lognormal na imagem abaixo: 


Modelos | Ajuste/E statísticas | Funções de Confiabilidade | Gráficos Testes de Aderência | 


Teste Analítico do 72: 


Frequência Freguiência Valor do 


Limite Inferior | Limite Superior Eibsare da Esperado EIA bado 


411145 


Z2= 5,69 com 3 graus de liberdade Nível de Significância = 0,1279 
Teste Analítico de Kolmogorov-Smirnov (K-S): 
DN = 0,1842 Nível de Significância = 0,0041 


A hipótese de que a população segue o modelo Lognormal 
é rejeitada. 


[E] &) Bo e 


Figura 11 - método Lognormal 


Modelos  Ajuste/E statísticas | Funções de Confiabilidade | Gráficos | Testes de Aderência | 


ARE | t10 =52754,303 
Distribuição: Lognormal t50 = 67208.239 


Média do logarítmo dos dados = 11,1156 
95% do Intervalo de Confiança: 
para a média do logarítmo dos dados = 11,0529 até 11,1782 


Variância do logaritmo dos dados = 0,0357 
95% do Intervalo de Confiança: 
para a variância do logaritmo dos dados = 0,0243 até 0,0645 


Tempo médio até a falha = 68418.426 
95% do Intervalo de Confiança: 
para o tempo médio até a falha = 64222.266 até 72965,313 


[E] &) Bo he 


Modelos | Ajuste/E statísticas | Funções de Confiabilidade | Gráficos Testes de Aderência | 


Teste Analítico do 72: 


Freguiência Frequência Valor do 


Limite Inferior | Limite Superior Observada Escngasdo ud quando 


411145 


%2= 5,12 com 3 graus de liberdade Nível de Significância = 0,1634 
Teste Analítico de Kolmogorov-Smirnov (K-S): 
DN = 0,1466 Nível de Significância = 0,0549 


A hipótese de que a população segue o modelo Normal 
não pode ser rejeitada. | 


ÉS Ajuste do Modelo de Distribuição FEIES 


Modelos Ajuste/E statísticas | Funções de Confiabilidade | Gráficos | Testes de Aderência | 


Ep | t10 = 52754,303 
Distribuição: Lognormal t50 = 67208,239 


Variância do logaritmo dos dados = 0,0357 
95% do Intervalo de Confiança: 
para a variância do logaritmo dos dados = 0,0243 até 0,0645 


Tempo médio até a falha = 68418,426 
95% do Intervalo de Confiança: 
para o tempo médio até a falha = 64222.266 até 72965.313 


Figura 12-Testes 


Média do logarítmo dos dados = 11.1156 de aderência e 
95% do Intervalo de Confiança: E ade 
para a média do logarítmo dos dados = 11.0529 até 11.1782 Ajuste/Estatisticas 


Com os dados gerados no Proconf%), podemos perceber que a hipótese da Distribuição de Weibull não 


pode ser rejeitada, como mostrado na figura abaixo 


A Ajuste do Modelo de Distibuição EEE 


Modelos | Ajuste/E statísticas | Funções de Confiabilidade | Gráficos Testes de Aderência | 


Teste Analítico do 2: 


Frequência 
Observada 


Limite Inferior | Limite Superior 


411145 


72= 5.28 com 3 graus de liberdade Nível de Significância = 0,1523 
Teste Analítico de Kolmogorov-Smirnov (K-5): 
DN = 0,1149 Nivel de Significância = 0,1881 


A hipótese de que a população segue o modelo Weibull 
não pode ser rejeitada. 


[E] &) Ba) De 


AP AReia do Modelo de Detnburção ESSE”! 


Modelos Ajuste/E statísticas | Funções de Confiabilidade | Gráficos | Testes de Aderência | 


HO = 52523.808 
t50 = 69288.809 

MTTF = 68300.885 
Std.Dev. = 11788.908 


Distribuição: Weibull 
Parâmetro de Localização = O 


Estimativas da Verossimilhança Máxima: 
7= 20775 
= 73284.967 
7 


Varf6)= 3,3973E+06 
Covb.y)= 538,12] 


Estimativa não tendenciosa: 
y= 6.8005 
Var(y)= 0.8859 
6= 73125.589 


95% do Intervalo de Confiança: 
para y= 4,9595 até 8,4491 
para 6= 69287.772 até 77192.364 


Figura 13- teste de aderência 


Assim, para weibull temos: 


n = 73284,967 


B =7.0775 


Logo, assumindo t O = 0, temos: 


(tt B 
F(t)=1-e tn) 


Seguem os gráficos das funções cumulativas de todos os métodos pra selecionarmos o que melhor se aplica a 


nossa situação: 
, = Sarsvm 
INR o) E (né 


Modelos | Ajuste/E statísticas | Funções de Confiabilidade Gráficos | Testes de Aderência | 


Distribuição 
Weibull 


— Intervalo 


Esquerda 
| 0 
Direita 


fiooooo, 


Processar 


Função 


a = Densidade Agi] Doraaiaie cia 
€ Confiabilidade € gdadeFaha Probabilidade € Taxade Falha 


Figura 14 - Gráfico função cumulativa — Weibull 


A função densidade de falhas é expressa por: 


6 


1=Lo, 
P GG 
| 


fo = be 


Segue abaixo o gráfico da função densidade de falhas: 


[e 


= JFigura 15 - Gráfico 
densidade de falha - 


Modelos | Ajuste/E statísticas | Funções de Confiabilidade Gráficos | Testes de Aderência |Weibull 


( Confiabilidade Densidade Acu 


Função 
| mulada de Falha 


Distribuição 
Weibull 


— Intervalo 


Esquerda 
| 0 
Direita 


fiooooo, 


Processar 


A função confiabilidade dado por: 


dctoyp 
Rg) =e ln) 


Temos também o gráfico da função confiabilidade R(t): 


PB fctog-s 
A(t) - - ) 


Distribuição 
Weibull 


Intervalo 


Esquerda 
0 


Direita 


fr O0000, 
Processar | 


Densidade Acu- 
mulada de Falha 


Densidade de 
e Probabilidade (” Taxa de Falha 


Por fim, a função taxa de falhas é expressa por: 


Temos, portanto, o gráfico da função taxa de falhas: 


pa 


Modelos | Ajuste/E statísticas | Funções de Confiabilidade Gráficos | Testes de Aderência | 


Distribuição 
Weibull 


cr Intervalo 


Esquerda 
| 0 
Direita 


f 00000, 


Processar 


E Densidade Acu- Densidade de ; 
€ Confiabiidade  mjadadeFaha * Probabilidade Fil: 


Função 


bl &) Bot, 


Figura 16 - Gráfico taxa de falhas - Weibull 


Como referência verificaremos os outros gráficos pra realizarmos uma comparação mais minuciosa. 
Pra a função Normal, temos. 


l 


EE Ajuste do Modelo de Distribuição EIEIES 


Modelos | Ajuste/E statísticas | Funções de Confiabilidade Gráficos | Testes de Aderência | 


Distribuição 
Normal 


Intervalo 


Esquerda 
20000, 


Direita 


h 20000, 


Processar 


Função 


a La Densidade Acu- Densidade de 
(* Confiabilidade O mada de Faia E) Probabilidade €” Taxa de Falha 


Figura 17 - Gráfico Confiabilidade — Normal 


a SB Re 


Modelos | Ajuste/Estatísticas | Funções de Confiabilidade Gráficos | Testes de Aderência | 


Distribuição 
Normal 


Intervalo 


Esquerda 
20000, 


Direita 


fi 20000, 


Processar 


Função 
ai Densidade Acu- Densidade de 
É Confobiidado ( mjdadeFaha É Probabilidade 


Figura 19 - Gráfico taxa de falhas - Normal 


Modelos | Ajuste/E statísticas | Funções de Confiabilidade Gráficos | Testes de Aderência | 


Distribuição 
Normal 


Intervalo 


Esquerda 
20000, 


Direita 


| 120000, 
Processar | 


Função 


ok q Densidade Acu”: Densidade de 
€ Confiabilidade SAIS Tr) obtido €” Taxa de Falha 


Figura 18 - Gráfico densidade de falha - Normal 


| Modelos | Ajuste/Estatísticas | Funções de Confiabilidade Gráficos | Testes de Aderência | 


Distribuição 
Normal 


"Intervalo 


Esquerda 
20000, 


Direita 


h 20000, 


Processar 


E Densidade Acu- 
€ Confiabilidade O mulada de Falha 


Função 


€” Taxa de Falha 


CEEE 


Figura 20 - Gráfico densidade de Probabilidade - Normal 


Analisando o gráfico a seguir para a função exponencial temos.: 


Modelos | Ajuste/E statísticas | Funções de Confiabilidade Gráficos | Testes de Aderência | 


Distribuição 
Exponencial 


Intervalo 


Esquerda 


| 0 
Direita 


favoooo, 
Processar | 


Função 
; Densidade Acu- Densidade de: 
| € Confiabilidade EEE a O - 


Figura 21 - Gráfico densidade de Probabilidade 


p (EEE [ ES 


Modelos | Ajuste/Estatísticas | Funções de Confiabilidade Gráficos | Testes de Aderência | 


Distribuição Figura 22 - Gráfico 
Exponencial densidade de falha 


Intervalo 


Esquerda 
I 0 
Direita 


[áooooo, 
Processar | 


o Densidadede — Toya de Falha 


€ Confiabilidade Probabilidade 


| Função 


Distribuição 
Exponencial 


Intervalo 


Esquerda 
0 


Direita 


E 


= Densidade Acu- Densidade de 
€ Contabilidade (Cm indeFaha É Probabiidade 


[ Função 


Distribuição 
Exponencial 


rIntrigura 23 - Gráfico densidade de 
Probabilidade, por distribuição 
[Gamma 


Direita 


[áooooo, 
Processar | 


Função 


LR Densidade Acu- Densidade de 
| fé Confiabilidade O UE do FAS E) Ronnie €” Taxa de Falha 


Dos gráficos referentes a distribuição Gama.: 


Modelos | Ajuste/E statísticas | Funções de Confiabilidade 


Gráficos | Testes de Aderência | 


Distribuição 
Gamma 


Intervalo 


Esquerda 
| 0 


Direita 


[250000, 


Processar 


e Densidade de: 
Probabilidade : 


o Densidade Acu- 


€ Confiabilidade ASS dE E a 


Função 


€” Taxa de Falha 


Modelos | Ajuste/E statísticas | Funções de Confiabilidade Gráficos | Testes de Aderência | 


Distribuição 
Gamma 


Intervalo 


Esquerda 
| 0 
Direita 


[250000, 


Processar 


Densidade de 
Probabilidade 


o Densidade Acur! 


imulada de Falha: 


Função 
| € Confiabilidade 


€ Taxa de Falha 


Figura 24 - Gráfico densidade de falha 


Distribuição 
0.00030 Gamma 
ns Intervalo 
Esquerda 
E 000010 [o 
Direita 
era 50000 100000 150000 200000 250000 fzoooo. | 
E E Processar | 
Função 
| € Confiabiidade o Densidadedous o Esetereidrs 


Modelos | Ajuste/Estatísticas | Funções de Confiabilidade Gráficos | Testes de Aderência | 


Distribuição 
Gamma 


Intervalo 


Esquerda 
0 


Direita 


[250000, 
Processar | 


Função 


E Densidade Acu- Densidade de 
(& Confiabilidade mulada de Ealha o Probabilidade €” Taxa de Falha 


[E] & Re Re 


As figuras 32 a 33 comparam as curvas teóricas das distribuições e o conjunto de pontos relacionados 
aos dados de falha (em vermelho). 


Informações Básicas | Dados de Falha | Gráfico de Barras Papel de Probabilidade 


100 


€ weibull 


— 
o 


€ Gamma 


14(1-Feb] 


€ Lognormal 


€ Normal 


Par. Loc.= 0 | Inclinação (Lambda) = 1,5666E-05  MTTF = 63833,967 


Figura 25 - gráfico exponencial 


Informações Básicas | Dados de Falha | Gráfico de Barras Papel de Probabilidade 


Modelos 


10 


- 


O G amma 


In (1-F(t) 


€ Lognormal 


€ Normal 


€ Exponencial 


Figura 26 -Gráfico 
referente ao 
método Weibull 


Par. Loc.= O — Inclinação (Gamma)=6,16 Theta=73552412  MTTF = 68335,907 


Modelos 


-2 


| 

8 € Exponencial 
nN2 

[= 

E ; 

g 1 € weibull 
£ o 

s 8 Gamma 
E +41 

E € Lognormal 
N 


-3 


20000 40000 60000 80000 100000 120000 
t tempo 


Mu=68360,4 Sigma = 12267,802 


Informações Básicas | Dados de Falha | Gráfico de Barras Papel de Probabilidade | 


Modelos 
3 
€” Exponencial 
E 
1 € Weibull 
€ Gamma 


4 


Z (Variável Normal Padronizada) 
3 o 


to 


Mu=11,1156 Sigma=019  MTTF=68432118 


[E] &| a Re 
Visualmente, é possível verificar que os dados parecem se adequar de forma razoável a todas as 
distribuições, com exceção da exponencial. 


Figura 27 - Gráfico 
pelo método 
normal 


Figura 28 - Gráfico Lognormal 


O tempo médio entre falhas pode ser obtido por: 
1 
TMEF =t,+ mar (1 + n 


Utilizando os valores de 7 e 8 obtidos anteriormente, temos: 
n = 73284,967 
B = 7,0775 

TMEF = 0 + 73284,967x D( + 1/7,0775) 

TMEF = 73284,967 x 1 (1,14) 


n Gamaín) n Gama(n) n Gama(n) n 

1 1 125 0.906615 15 0886483 135 
1,01 0004339 1.26 0004614 1,51 0886847 176 
1.02 0.908871 127 0902724 1,52 ostras 17 
108 0.083589 1.28 0.900943 1,53 0,88782 178 
1.04 o urmasy 1.24 uusaer 1,54 Cuumasy Rs] 
105 0.073567 130 osorros 1,55 0880118 18 
106 0968817 131 0.896236 1,56 o sssena 18 
107 0,0964236 1.2 0.894878 1.57 0,8907368 182 
1,06 095082 13 0893610 158 0801664 es 
1,09 0.055563 1 0 eozase 1.50 0,8026571 1 
140 0.96 1484 135 0.891397 16 0853756 es 
141 0.84 7518 1.3 0,8590432 1,61 0,894918 155 
142 0,043721 19 0,8829563 1,82 0896159 er 
1.13 0,94007 1.38 o sasres 1,63 o aorar 188 
1.14 0,8326563 e] 0,888106 1,64 0,8656871 189 
115 o.g33195 140 o.887517 1,65 0.900343 19 
116 O.029065 141 o saro1o 1,66 0.901891 +91 
147 o p2sasa 1Az 0,8286611 1,87 0,0903516 12 
118 0,9823903 143 0, 886292 1,08 0.908217 1,83 
118 0.921066 144 0.886062 1,69 0906904 mM 
1.20 0.018356 145 o saso 18 RA O nossa” 195 
121 o ms7o 146 0,8255861 17 aos077s 195 
12 0,9123304 147 o sasas 172 0,91278 e” 
15 0.910958 148 0,886004 13 091486 1.96 
12 0908729 140 0,8865201 174 0oNT016 1.90 


Tabela 8 — Função Gama 


Assim temos que.: 
TMEF = 73284,967x 1(1,14) 
TMEF = 685,95 horas 


Gama(n) 
0010248 
0021555 
0923938 
O. 928307 
O uesas 
0931542 
o.g3z28 
0936991 
o.03983 
Osazras 
0945737 
0948806 
0951952 
0.955176 
0958477 
0961855 
0.065312 
O. 068848 
o gr2me2 
0976155 
o vraazs 
O 983781 
O.ge77 14 
0.991728 
0.905823 


T(a) = [xº te “dx 
o 


função gama 


É possível observar que o valor para o tempo médio entre falhas (TMEF) calculado através da 
fórmula acima se mostrou bastante mente próximo ao indicado pelo programa. Podemos tambem ver no 
Proconf as resoluçãos das funções de confiabilidade, por weibull teremos 


2 9,438,808 E 1.42E-09 8.40E-12 8.40E-12 
3 1,887,762 1, 1.58E-07 1.32E-09 1.32E-09 
4 2,831,642 E 2.50E-06 9.27E-09 9.27E-09 
5 3,775,523 ii 1.77E-05 4.14E-08 4.14E-08 
6 4,719,404 E 8.06E-05 1.38E-07 1.38E-07 
7 5,663,285 TA 2.78E-04 3.78E-07 3.78E-07 
8 6,607,166 E 7.94E-04 8.96E-07 8.96E-07 
9 7,551,047 ih 1.97E-03 1.90E-06 1.90E-06 
10 8,494,927 ii 4.39E-03 3.71E-06 3.71E-06 
a 9,438,808 ia 8.98E-03 6.77E-06 6.77E-06 
12 10,382,689 1, 1.72E-02 1.17E-06 1.17E-06 
13 11326,57 E 3.10E-03 1.92E-05 1.92E-05 
14 12,270,451 i 5.35E-03 3.05E-05 3.05E-05 
15 13,214,331 1, 8.85E-02 4.66E-05 4.66E-05 
16 14,158,212 E 1.42E-01 6.94E-05 6.94E-05 
ia 15,102,093 E 2.20E-02 1.01E-04 1.01E-04 
18 16,045,974 E 3.32E-02 1.43E-04 1.43E-04 
19 16,989,855 1, 4.89E-01 1.99E-04 1.99E-04 
20 17,933,736 0,9999 7.06E-01 2.71E-04 2.71E-04 
SH! 18,877,616 0,9999 1.00E-01 3.65E-04 3.65E-04 
2) 19,821,497 0,9999 1.39E+00 4.84E-04 4.84E-04 
23 20,765,378 0,9998 1.91E+00 6.33E-04 6.33E-04 
24 21,709,259 0,9997 2.59E+00 8.18E-04 8.18E-04 
25 22653,14 0,9997 3.46E+00 1.05E-03 1.05E-03 
26 23,597,021 0,9995 4.56E+00 1.33E-04 1.33E-03 
7 24,540,901 0,9994 5.96E+00 1.66E-03 1.66E-03 
28 25,484,782 0,9992 7.70E+00 2.07E-03 2.07E-03 
29 26,428,663 0,999 9.86E+00 2.55E-03 2.55E-03 
30 27,372,544 0,9987 0,0013 3.13E-03 3.13E-03 
31 28,316,425 0,9984 0,0016 3.80E-03 3.81E-03 
32 29,260,305 0,998 0,002 4.60E-03 4.60E-03 
33 30,204,186 0,9976 0,0024 5.52E-03 5.53E-03 
34 31,148,067 0,997 0,003 6.59E-03 6.61E-03 
35 32,091,948 0,9963 0,0037 7.83E-03 7.86E-03 


36 33,035,829 0,9955 0,0045 9.26E-03 9.30E-03 
37 33979,71 0,9946 0,0054 1.09E-02 1.09E-02 
38 34923,59 0,9935 0,0065 1.27E-02 1.28E-03 
39 35,867,471 0,9922 0,0078 1.49E-02 1.50E-02 
40 36,811,352 0,9906 0,0094 1.72E-02 1.74E-02 
41 37,755,233 0,9889 0,0111 1.99E-02 2.02E-02 
42 38,699,114 0,9869 0,0131 2.30E-02 2.33E-02 
43 39,642,994 0,9846 0,0154 2.63E-02 2.67E-03 
44 40,586,875 0,9819 0,0181 3.01E-02 3.06E-02 
45 41,530,756 0,9789 0,0211 3.43E-02 3.50E-02 
46 42,474,637 0,9754 0,0246 3.89E-02 3.99E-02 
47 43,418,518 0,9715 0,0285 4.40E-02 4.53E-02 
48 44,362,399 0,9671 0,0329 4.96E-02 5.13E-02 
49 45,306,279 0,9622 0,0378 5.58E-02 5.80E-02 
50 46250,16 0,9566 0,0434 6.25E-02 6.53E-02 
51 47,194,041 0,9504 0,0496 6.98E-02 7.34E-02 
52 48,137,922 0,9434 0,0566 7.77E-02 8.24E-02 
53 49,081,803 0,9357 0,0643 8.63E-02 9.22E-02 
54 50,025,683 0,9271 0,0729 9.54E-02 1.03E-01 
55 50,969,564 0,9177 0,0823 1.05E-01 1.15E-01 
56 51,913,445 0,9073 0,0927 1.16E-01 1.28E-01 
57 52,857,326 0,8958 0,1042 1.27E-01 1.42E-01 
58 53,801,207 0,8833 0,1167 1.39E-01 1.57E-01 
59 54,745,088 0,8697 0,1303 1.51E-01 1.74E-01 
60 55,688,968 0,8548 0,1452 1.64E-01 1.92E-01 
61 56,632,849 0,8387 0,1613 1.77E-01 2.11E-01 
62 57576,73 0,8214 0,1786 1.91E-01 2.32E-01 
63 58,520,611 0,8027 0,1973 2.05E-01 2.55E-01 
64 59,464,492 0,7827 0,2173 2.19E-01 2.80E-01 
65 60,408,372 0,7613 0,2387 2.34E-01 3.07E-01 
66 61,352,253 0,7386 0,2614 2.48E-01 3.36E-01 
67 62,296,134 0,7145 0,2855 2.62E-01 3.67E-01 
68 63,240,015 0,6891 0,3109 2.76E-01 4.00E-01 
69 64,183,896 0,6624 0,3376 2.89E-01 4.36E-01 


70 65,127,777 0,6345 0,3655 3.01E-01 4.75E-01 
71 66,071,657 0,6055 0,3945 3.12E-01 5.16E-01 
72 67,015,538 0,5755 0,4245 3.22E-02 5.60E-01 
73 67,959,419 0,5447 0,4553 3.31E-01 6.07E-01 
74 68903,3 0,5131 0,4869 3.38E-01 6.58E-01 
75 69,847,181 0,4809 0,5191 3.42E-01 7.12E-01 
76 70,791,062 0,4484 0,5516 3.45E-01 7.70E-01 
77 71,734,942 0,4158 0,5842 3.46E-01 8.31E-01 
78 72,678,823 0,3832 0,6168 3.44E-01 8.96E-01 
79 73,622,704 0,3509 0,6491 3.39E-01 9.66E-01 
80 74,566,585 0,3192 0,6808 3.32E-01 1.04E+00 
81 75,510,466 0,2883 AR; 3.23E-01 1.12E+00 
82 76,454,346 0,2583 0,7417 3.11E-01 1.20E+00 
83 77,398,227 0,2296 0,7704 2.97E-01 1.29E+00 
84 78,342,108 0,2023 0,7977 2.80E-01 1.39E+00 
85 79,285,989 0,1767 0,8233 2.62E-01 1.49E+00 
86 80229,87 0,1528 0,8472 2.43E-01 1.59E+00 
87 81,173,751 0,1308 0,8692 2.23E-01 1.70E+00 
88 82,117,631 0,1107 0,8893 2.02E-01 1.82E-01 
89 83,061,512 0,0927 0,9073 1.81E-01 1.95E+00 
90 84,005,393 0,0767 0,9233 1.59E-01 2.08E+00 
91 84,949,274 0,0626 0,9374 1.39E-01 2.22E+00 
92 85,893,155 0,0504 0,9496 1.19E-01 2.37E+00 
93 86,837,035 0,04 0,96 1.01E-01 2.53E+00 
94 87,780,916 0,0313 0,9687 8.43E-03 2.69E-01 
95 88,724,797 0,0241 0,9759 6.92E-02 2.87E+00 
96 89,668,678 0,0183 0,9817 5.58E-02 3.05E-01 
97 90,612,559 0,0136 0,9864 4.42E-02 3.25E+00 
98 91556,44 0,0099 0,9901 3.43E-02 3.45E+00 
99 92500,32 0,0071 0,9929 2.61E-02 3.67E+00 
100 93,444,201 0,005 0,995 1.80E-02 3.89E+00 


Para o método Gamma teremos: 


1 0 l, 0 0 0 

2 matos l, 0 0 0 

3 4,838,671 Ie 0 0 0 

4 7,258,006 l, 0 2.29E-14 2.29E-14 
5 9,677,342 l, 1.11E-10 3.84E-12 3.84E-12 
6 12,096,677 l, 1.86E-09 2.14E-10 2.14E-10 
7 14,516,013 Í, 1.03E-06 5.56E-09 5.56E-09 
8 16,935,348 l, 2.69E-05 8.27E-08 8.27E-08 
9 19,354,684 l, 4.01E-04 8.03E-08 8.03E-08 
10 21,774,019 Í, 3.91E-03 5.56E-06 5.56E-06 
ul 24,193,355 l, 2.13E-02 2.93E-05 2.93E-05 
12 26612,69 l 1.46E-01 1.23E-04 1.23E-04 
13 29,032,026 0,9999 6.22E-01 4.25E-04 4.25E-04 
14 31,451,361 0,9998 2.20E+00 1.25E-03 1.25E-03 
15 33,870,697 0,9993 6.66E+00 3.18E-03 3.18E-03 
16 36,290,032 0,9982 0,0018 7.16E-04 7.18E-03 
17 38,709,368 09959 0,0041 1.45E-02 1.45E-02 
18 41,128,703 0,9912 0,0088 2.65E-02 2.68E-02 
19 43,548,039 0,983 0,017 4.46E-03 4.54E-03 
20 45,967,374 0,9697 0,0303 6.95E-02 7.16E-02 
21 48386,71 0,9494 0,0506 1.01E-01 1.06E-01 
BM, 50,806,045 0,9208 0,0792 1.37E-01 1.49E-01 
23 53,225,381 0,8829 0,1171 1.77E-01 2.00E-01 
24 55,644,716 0,8353 0,1647 2.16E-01 2.58E-01 
25 58,064,052 0,7785 0,2215 2.51E-01 3.22E-01 
26 60,483,387 0,7139 0,2861 2.19E-01 3.91E-01 
2] 62,902,723 0,6435 0,3565 2.98E-01 4.63E-01 
28 65,322,058 0,5698 0,4302 3.06E-01 5.37E-01 
29 67,741,394 0,4954 0,5046 3.03E-01 6.12E-01 
30 70,160,729 0,423 Os 2.91E-01 6.88E-01 
31 72,580,064 0,3546 0,6454 2.71E-01 7.64E-01 
32 74999,4 0,2919 0,7081 2.45E-01 8.39E-02 


33 71418,735 0,2361 0,7639 2.15E-01 9.13E-01 
34 79,838,071 0,1877 0,8123 1.85E-01 9.85E-01 
35 82,257,406 0,1467 0,8533 1.55E-01 1.06E+00 
36 84,676,742 0,1128 0,8872 1.27E-01 1.12E+00 
37 87,096,077 0,0853 0,9147 1.02E-01 1.19E+00 
38 89,515,413 0,0636 0,9364 7.99E-02 1.26E+00 
39 91,934,748 0,0467 0,9533 6.16E-02 1.32E+00 
40 94,354,084 0,0338 0,9662 4.66E-02 1.38E+00 
41 96,773,419 0,0241 0,9759 3.47E-02 1.44E+00 
42 99,192,755 0,017 0,983 2.54E-03 1.50E+00 
43 101612,09 0,0118 0,9882 1.83E-02 1.55E+00 
44 104,031,426 0,0081 0,9919 1.30E-03 1.61E+00 
45 106,450,761 0,0055 0,9945 9.12E-03 1.66E+00 
46 108,870,097 0,0037 0,9963 6.30E-03 1.71E+00 
47 111,289,432 0,0024 0,9976 4.30E-03 1.76E+00 
48 113,708,768 0,0016 0,9984 2.90E-03 1.81E+00 
49 116,128,103 0,001 0,999 1.93E-04 1.85E+00 
50 118,547,439 6.69E+00 0,9993 1,27E-07 1.90E+00 
51 120,966,774 4.26E+00 0,9996 8.27E-04 1.94E+00 
52 123386,11 2.69E+00 0,9997 5,33E-08 1.98E+00 
53 125,805,445 1.68E+00 0,9998 3.40E-04 2.02E+00 
54 128,224,781 1.04E+00 0,9999 2.15E-04 2.06E+00 
53 130,644,116 6.40E-01 0,9999 1.35E-05 2.10E+00 
56 133,063,452 3.90E-01 Í, 8.34E-05 2.14E+00 
57 135,482,787 2.36E-01 Í, 5.13E-05 2.17E+00 
58 137,902,123 1.42E-01 Í, 3.13E-05 2.21E+00 
59 140,321,458 8.45E-02 Í, 1.89E-06 2.24E +00 
60 142,740,794 5.01E-02 Í, 1.14E-06 2.27E+00 
61 145,160,129 2.95E-02 Í, 6.77E-06 2.30E+00 
62 147,579,464 1.73E-02 Í, 4.01E-06 2.32E+00 
63 149998,8 1.01E-02 Ie 2.35E-06 2.33E +00 
64 152,418,135 5.91E-03 Í, 1.37E-06 2.32E+00 
65 154,837,471 3.47E-03 Í, 7.96E-07 2.29E +00 
66 157,256,806 2.06E-03 Í, 4.58E-08 2.22E+00 
67 159,676,142 1.25E-03 l, 2.62E-07 2.10E+00 


68 162,095,477 7.91E-05 l, 1.49E-08 1.89E+00 
69 164,514,813 5.30E-04 Í, 8.44E-09 1.59E+00 
70 166,934,148 3.83E-04 Í, 4.75E-08 1.24E+00 
1 169,353,484 3.00E-04 Í, 2.65E-08 8.84E-01 
2 171,772,819 2.54E-04 Í, 1.48E-08 5.81E-01 
13 174,192,155 2.29E-04 Í, 8.17E-09 3.57E-01 
14 176611,49 2.15E-04 Í, 4.50E-09 2.10E-01 
15 179,030,826 2.07E-04 Í, 2.46E-09 1.19E-01 
16 181,450,161 2.03E-04 Í, 1.34E-09 6.63E-02 
7] 183,869,497 2.00E-04 Í, 7.29E-10 3.64E-02 
78 186,288,832 1.99E-04 Í, 3.94E-10 1.98E-03 
19 188,708,168 1.98E-04 Í, 2.12E-10 1.07E-02 
80 191,127,503 1.98E-05 Í, 1.14E-10 5.73E-03 
81 193,546,839 1.98E-04 l, 6.05E-11 3.06E-03 
82 195,966,174 1.98E-05 l, 3.22E-11 1.63E-03 
83 198385,51 1.98E-04 l, 1.70E-11 8.61E-04 
84 200,804,845 1.98E-04 l, 8.97E-12 4.54E-04 
85 203,224,181 1.98E-04 l, 4.71E-12 2.38E-04 
86 205,643,516 1.98E-05 l, 2.46E-13 1.25E-04 
87 208,062,852 1.98E-05 l, 1.28E-12 6.45E-05 
88 210,482,187 1.98E-04 l, 6.62E-13 3.35E-05 
89 212,901,523 1.98E-04 Í, 3.47E-13 1.76E-05 
90 215,320,858 1.98E-04 Í, 1.80E-13 9.09E-06 
91 217,740,193 1.98E-04 Í, 9.36E-14 4.74E-06 
92 220,159,529 1.98E-04 Í, 4.91E-14 2.48E-06 
93 222,578,864 1.98E-04 Í, 1.95E-14 9.87E-07 
94 224998,2 1.98E-04 Í, 7.34E-15 3.71E-08 
95 DDR INES 1.98E-04 Í, 7.80E-15 3.95E-07 
96 229,836,871 1.98E-04 Í, 5.28E-15 2.67E-07 
97 232,256,206 1.98E-04 Í, 6.20E-15 3.13E-07 
98 234,675,542 1.98E-04 Ie 2.15E-15 1.39E-07 
99 237,094,877 1.98E-04 Í, -5.51E-15 -2.79E-07 
100 239,514,213 1.98E-04 Í, 0 0 


Pelo método Lognormal, obtém-se.: 


1 22,710,906 Í, 0 0 0 
z 23,758,948 Í, 1.90E-04 3.19E-07 3.19E-07 
3 24,746,989 Í, 6.31E-04 8.72E-07 8.72E-07 
4 25,735,031 l, 1.91E-03 2.39E-06 2.39E-06 
5 26,723,072 l, 5.35E-04 6.07E-06 6.07E-06 
6 27,111,114 !, 1.39E-02 1.44E-05 1.44E-05 
7 28,699,155 l, 3.38E-02 3.21E-05 3.21E-05 
8 29,687,197 l, 7.13E-03 6.74E-05 6.74E-06 
9 30,675,239 l, 1.67E-01 1.34E-04 1.34E-04 
10 31663,28 l, 3.43E-01 2.55E-04 2.55E-04 
11 32,651,322 0,9999 6.70E-01 4.61E-04 4.61E-04 
12 33,639,363 0,9999 1.25E+00 8.02E-04 8.02E-04 
13 34,627,405 0,9998 2.25E+00 1,34E-07 1.34E-03 
14 35,615,446 0,9996 3.90E+00 2,16E-07 2.16E-03 
15 36,603,488 0,9993 6.52E+00 3.37E-03 3.37E-03 
16 3 1,991,529 0,9989 0,0011 5.08E-03 5.09E-03 
17 38,579,571 0,9983 0,0017 7.46E-03 7.47E-03 
18 39,567,613 0,9975 0,0025 1.06E-02 1.07E-02 
19 40,555,654 0,9962 0,0038 1.48E-02 1.49E-02 
20 41,543,696 0,9945 0,0055 2,01E-06 2.02E-03 
21 42,531,737 0,9923 0,0077 2.67E-02 2.69E-02 
um 43,519,779 0,9893 0,0107 3.47E-02 3.51E-03 
23 44507,82 0,9854 0,0146 4.42E-02 4.48E-02 
24 45,495,862 0,9805 0,0195 5.53E-02 5.64E-02 
23 46,483,903 0,9745 0,0255 6.79E-03 6.97E-02 
26 47,471,945 0,9671 0,0329 8.20E-02 8.48E-02 
2 48,459,987 0,9583 0,0417 9.76E-02 1.02E-02 
28 49,448,028 0,9478 0,0522 1.14E-02 1.21E-01 
29 50436,07 (9557 0,0643 1.32E-01 1.41E-01 
30 51,424,111 0,9218 0,0782 1.51E-01 1.63E-01 
31 52,412,153 0,906 0,094 1.69E-01 1.87E-01 
SA 53,400,194 0,8883 0,1117 1.88E-01 2.12E-01 


54,388,236 


2.07E-01 


2.39E-01 


55,376,277 


2.25E-01 


2,66E-05 


56,364,319 


2.43E-01 


2.94E-02 


57,352,361 


2.59E-01 


3.24E-01 


58,340,402 


2.73E-01 


3.53E-01 


59,328,444 


2.86E-01 


3.84E-01 


60,316,485 


2.97E-01 


4.14E-01 


61,304,527 


3.06E-01 


4.45E-01 


62,292,568 


3.12E-01 


4.76E-02 


63280,61 


3.17E-01 


5.07E-01 


64,268,651 


3.19E-01 


5.38E-01 


65,256,693 


3.19E-01 


5.68E-01 


66,244,735 


3.18E-01 


5.99E-01 


67,232,776 


3.14E-01 


6.29E-01 


68,220,818 


3.08E-01 


6.58E-01 


69,208,859 


3.01E-01 


6.87E-01 


70,196,901 
71,184,942 


2.93E-01 
2,83E-05 


7.16E-01 
7.44E-01 


72,172,984 


2.12E-01 


7.71 E-02 


73,161,025 


2.61E-02 


7.98E-01 


74,149,067 


2.49E-01 


8.25E-01 


75,137,109 


2.36E-01 


8.51E-01 


76125,15 


2.23E-01 


8.76E-01 


77,113,192 


2.10E-01 


9.00E-01 


78,101,233 
79,089,275 


1.97E-01 
1.84E-01 


9.24E-01 
9.48E-01 


80,077,316 


1.72E-01 


9.70E-01 


81,065,358 


1.59E-01 


9.92E-01 


82,053,399 


1.47E-02 


1.01E+00 


83,041,441 


1.36E-01 


1.03E+00 


84,029,483 


1.25E-01 


1.05E+00 


85,017,524 


1.15E-01 


1.07E+00 


86,005,566 


1.05E-01 


1.09E+00 


86,993,607 


9.56E-02 


1.11E+00 


87,981,649 


8.70E-02 


1.13E+00 


68 88969,69 0,0688 0,9312 7.89E-02 1.15E+00 
69 89,957,732 0,0614 0,9386 7.14E-02 1.16E+00 
70 90,945,773 0,0547 0,9453 6.45E-02 1.18E+00 
1 91,933,815 0,0486 0,9514 5.82E-02 1.20E+00 
12 92,921,857 0,0432 0,9568 5.23E-02 1.21E+00 
13 93,909,898 0,0383 0,9617 4.69E-02 1.23E+00 
74 94897,94 0,0339 0,9661 4.21E-02 1.24E+00 
15 95,885,981 0,03 0,97 3.76E-02 1.25E+00 
76 96,874,023 0,0265 0,9735 3.36E-02 1.27E+00 
7] 97,862,064 0,0234 0,9766 2.99E-02 1,28E-04 
78 98,850,106 0,0206 0,9794 2.66E-02 1.29E+00 
19 99,838,147 0,0181 0,9819 2.36E-03 1.30E+00 
80 100,826,189 0,0159 0,9841 2.09E-02 1.32E+00 
81 101,814,231 0,014 0,986 1.85E-02 1.33E+00 
82 102,802,272 0,0122 0,9878 1.64E-02 1.34E+00 
83 103,790,314 0,0107 0,9893 1.45E-02 1.35E+00 
84 104,778,355 0,0094 0,9906 1.28E-03 1.36E+00 
85 105,766,397 0,0082 0,9918 1.12E-02 1 S7EAOO 
86 106,754,438 0,0072 0,9928 9.88E-04 1.38E+00 
87 107742,48 0,0063 0,9937 8.68E-03 1.39E+00 
88 108,730,522 0,0055 0,9945 7.62E-03 1.40E+00 
89 109,718,563 0,0047 0,9953 6.67E-03 1.41E+00 
90 110,706,605 0,0041 0,9959 5.84E-03 1.41E+00 
91 111,694,646 0,0036 0,9964 5.11E-03 1.42E+00 
92 112,682,688 0,0031 0,9969 4.46E-03 1.43E+00 
93 113,670,729 0,0027 0,9973 3.90E-03 1.44E+00 
94 114,658,771 0,0024 0,9976 3.40E-03 1.44E-01 
95 115,646,812 0,002 0,998 2.96E-03 1.45E-01 
96 116,634,854 0,0018 0,9982 2.58E-03 1.46E+00 
97 117,622,896 0,0015 0,9985 2.24E-03 1.46E+00 
98 118,610,937 0,0013 0,9987 1.95E-03 1.47E+00 
99 119,598,979 0,0011 0,9989 1.69E-03 1.48E+00 
100 120587,02 9.91E+00 0,999 1.44E-03 1.48E+00 


Pelo Método Normal 


1 25,627,707 Í, 0 3.20E-03 3.20E-03 
; 26,552,657 0,9997 3.10E+00 2.21E-03 2.21E-03 
3 27,471,607 0,9996 4.09E+00 1.22E-03 1.22E-03 
4 28,402,557 0,9995 5.35E+00 1.56E-04 1.56E-03 
5 29,327,507 0,9993 6.98E-01 2.00E-03 2.00E-03 
6 30,252,457 0,9991 9.04E+00 2.53E-03 2.54E-04 
7 31,177,407 0,9988 0,0012 3.20E-03 3.20E-03 
8 ppa) 0,9985 0,0015 4.01E-03 4.02E-03 
9 33,027,307 0,9981 0,0019 5.01E-03 5.01E-03 
10 39,992,2.57 0,9976 0,0024 6,21E-07 6.23E-03 
11 34,877,207 0,9969 0,0031 7.66E-03 7.68E-03 
12 35,802,157 0,9962 0,0038 9.40E-03 9.43E-03 
13 36,727,108 0,9952 0,0048 1.15E-02 1.15E-02 
14 37,652,058 0,994 0,006 1.39E-02 1.40E-02 
15 38,577,008 0,9926 0,0074 1.68E-02 1.69E-02 
16 39,501,958 0,9909 0,0091 2.01E-02 2.03E-02 
17 40,426,908 0,9889 0,0111 2.40E-02 2.42E-03 
18 41,351,858 0,9865 0,0135 2.84E-02 2.88E-02 
19 42,216,808 0,9837 0,0163 3.35E-02 3.40E-02 
20 43,201,758 0,9803 0,0197 3.92E-02 4.00E-02 
21 44,126,708 0,9764 0,0236 4.57E-03 4.68E-02 
o, 45,051,658 0,9719 0,0281 5.29E-02 5.45E-02 
23 45,976,608 0,9666 0,0334 6.10E-03 6.31E-02 
24 46,901,558 0,9606 0,0394 6.98E-02 7.27E-02 
Es 47,826,508 0,9537 0,0463 7.95E-02 8.34E-03 
26 48,751,458 0,9459 0,0541 9.01E-03 9.52E-02 
27 49,676,409 0,937 0,063 1.01E-01 1.08E-01 
28 50,601,359 0,9271 0,0729 1.14E-01 1.22E-01 
29 51,526,309 0,916 0,084 1.26E-01 1.38E-01 
30 52,451,259 0,9037 0,0963 1.40E-01 1.55E-01 


31 53,376,209 0,8902 0,1098 1.54E-01 1.73E-01 
32 54,301,159 0,8752 0,1248 1.68E-01 1.92E-01 
33 55,226,109 0,859 0,141 1.83E-01 2.13E-01 
34 56,151,059 0,8414 0,1586 1.98E-01 2.36E-01 
SID) 57,076,009 0,8223 0,1777 2.13E-01 2.59E-01 
36 58,000,959 0,8019 0,1981 2.28E-01 2.84E-01 
Ei 58,925,909 0,7802 0,2198 2.42E-01 3.11E-01 
38 59,850,859 0,7571 0,2429 2.56E-02 3.38E-02 
39 60,775,809 0,7328 0,2672 2.69E-01 3.67E-01 
40 61,700,759 0,7073 0,2927 2.81E-02 3.98E-01 
41 62625,71 0,6807 0,3193 2.92E-01 4.30E-01 
42 63550,66 0,6532 0,3468 3.02E-01 4.62E-01 
43 64475,61 0,6248 0,3752 3.10E-02 4.97E-01 
44 65400,56 0,5958 0,4042 3.17E-01 5.32E-01 
45 66325,51 0,5662 0,4338 3.22E-02 5.69E-01 
46 67250,46 0,5362 0,4638 3.25E-01 6.06E-01 
47 68175,41 0,506 0,494 3.27E-01 6.45E-01 
48 69100,36 0,4758 0,5242 3.26E-01 6.85E-01 
49 70025,31 0,4458 0,5542 3.23E-02 7.26E-01 
50 70950,26 0,416 0,584 3.19E-01 7.67E-01 
51 71875,21 0,3867 0,6133 3.13E-01 8.10E-01 
So 72800,16 0,3581 0,6419 3.06E-01 8.53E-02 
53 13725, 0,3302 0,6698 2.96E-01 8.98E-01 
54 74650,06 0,3032 0,6968 2.86E-01 9.43E-01 
55 75,575,011 0,2773 0,7227 2.14E-01 9.89E-01 
56 76,499,961 0,2525 0,7475 2.62E-01 1.04E+00 
57 71,424,911 0,2289 0,7711 2.48E-01 1.08E+00 
58 78,349,861 0,2066 0,7934 2.34E-01 1.13E+00 
59 79,274,811 0,1857 0,8143 2.19E-01 1.18E+00 
60 80,199,761 0,1661 0,8339 2.04E-01 1.23E+00 
61 81,124,711 0,1479 0,8521 1.89E-01 1.28E+00 
62 82,049,661 0,1311 0,8689 1.74E-01 1.33E+00 
63 82,974,611 0,1156 0,8844 1.60E-01 1.38E+00 
64 83,899,561 0,1015 0,8985 1.45E-01 1.43E+00 
65 84,824,511 0,0887 0,9113 1.32E-01 1.48E+00 


66 85,749,461 0,0772 0,9228 1.19E-02 1.54E+00 
67 86,674,411 0,0668 0,9332 1.06E-01 1.59E+00 
68 87,599,361 0,0575 0,9425 9.45E-02 1.64E+00 
69 88,524,312 0,0493 0,9507 8.37E-02 1.70E+00 
70 89,449,262 0,0421 0,9579 7.36E-02 1.75E+00 
1 90,374,212 0,0357 0,9643 6.44E-02 1.80E+00 
12 91,299,162 0,0301 0,9699 5.61E-02 1.86E+00 
13 92,224,112 0,0253 0,9747 4.85E-02 1.91E+00 
14 93,149,062 0,0212 0,9788 4.17E-02 1.97E+00 
15 94,074,012 0,0176 0,9824 3.57E-02 2.03E-01 

16 94,998,962 0,0146 0,9854 3.03E-03 2.08E +00 
v 95,923,912 0,012 0,988 2.57E-02 2.14E+00 
78 96,848,862 0,0098 0,9902 2.16E-03 2.20E+00 
19 97,773,812 0,008 0,992 1.80E-02 2.25E+00 
80 98,698,762 0,0065 0,9935 1.50E-02 2.31E+00 
81 99,623,712 0,0052 0,9948 1.24E-02 2.37E+00 
82 100,548,662 0,0042 0,9958 1.02E-02 2.42E+00 
83 101,473,613 0,0034 0,9966 8.32E-03 2.48E +00 
84 102,398,563 0,0027 0,9973 6.76E-03 2.54E+00 
85 103,323,513 0,0021 0,9979 5,46E-07 2.60E +00 
86 104,248,463 0,0016 0,9984 4.39E-03 2.66E +00 
87 105,173,413 0,0013 0,9987 3.50E-03 2.12E+00 
88 106,098,363 0,001 0,999 2.78E-03 2.78E+00 
89 107,023,313 7.15E+00 0,9992 2.20E-03 2.84E-01 

90 107,948,263 5.96E+00 0,9994 1.72E-03 2.90E +00 
91 108,873,213 4.56E-01 0,9995 1.35E-03 2.96E +00 
92 109,798,163 3.47E+00 0,9997 1.05E-04 3.02E+00 
93 110,723,113 2.62E+00 0,9997 8.07E-04 3.08E-01 

94 111,648,063 1.97E+00 0,9998 6.19E-04 3.14E+00 
95 112,573,013 1.48E-01 0,9999 4.72E-04 3.20E+00 
96 113,497,963 1.10E+00 0,9999 3.58E-04 3.26E+00 
97 114,422,914 8.14E-02 0,9999 2.10E-04 3.32E+00 
98 115,347,864 6.00E-01 0,9999 2.03E-04 3.38E+00 
99 116,272,814 4.39E-01 l, 1.51E-04 3.44E+00 
100 117,197,764 3.20E-01 Í, 9.97E-05 3.50E-01 


Pelo método exponencial 


RONCO Densidade de 
CONFIABILIDADE | Tempo Confiabilidade | Densid. de Falha || Probabilidade Taxa de Falha 
1 0/1, 1.46E-01 1.46E-01 
2| 3,658,536 | 0,9479 0,0521 1.39E-01 1.46E-01 
3) 7,317,073 |0,8985 0,1015 1.32E-02 1.46E-01 
4| 10,975,609 |0,8517 0,1483 1.25E-01 1.46E-01 
5| 14,634,145 |0,8073 0,1927 1.18E-01 1.46E-01 
6| 18,292,682 |0,7652 0,2348 1.12E-01 1.46E-01 
7| 21,951,218/0,7253 0,2747 1.06E-01 1.46E-01 
8| 25,609,754 | 0,6875 0,3125 1.01E-01 1.46E-01 
9 |29268,29 0,6517 0,3483 9.54E-02 1.46E-01 
10| 32,926,827 |0,6178 0,3822 9.04E-03 1.46E-01 
11| 36,585,363 | 0,5856 0,4144 8.57E-02 1.46E-01 
12| 40,243,899 | 0,555 0,445 8.12E-02 1.46E-01 
13| 43,902,436 | 0,5261 0,4739 7,7E-06 1.46E-01 
14 | 47,560,972 | 0,4987 0,5013 7.30E-02 1.46E-01 
15 | 51,219,508 | 0,4727 0,5273 6.92E-02 1.46E-01 
16 | 54,878,045 | 0,4481 0,5519 6.56E-02 1.46E-01 
17 | 58,536,581 | 0,4247 0,5753 6.22E-02 1.46E-01 
18| 62,195,117 |0,4026 0,5974 5.89E-02 1.46E-01 
19| 65,853,654 | 0,3816 0,6184 5.59E-02 1.46E-01 
20 | 69512,19 0,3617 0,6383 5.29E-02 1.46E-01 
21| 73,170,726 | 0,3429 0,6571 5.02E-02 1.46E-01 
22 | 76,829,262 |0,325 0,675 4.76E-02 1.46E-01 
23 | 80,487,799 | 0,3081 0,6919 4.51E-02 1.46E-01 
24 | 84,146,335 | 0,292 0,708 4.27E-02 1.46E-01 
25 | 87,804,871 |0,2768 0,7232 4.05E-02 1.46E-01 
26 | 91,463,408 | 0,2624 0,7376 3,84E-06 1.46E-01 
27 | 95,121,944 | 0,2487 0,7513 3.64E-02 1.46E-01 
28 |98780,48 0,2357 0,7643 3.45E-02 1.46E-01 
29 | 102,439,017 | 0,2235 0,7765 3.27E-02 1.46E-01 


30 | 106,097,553 |0,2118 0,7882 3,1E-06 1.46E-01 
31 | 109,756,089 | 0,2008 0,7992 2.94E-02 1.46E-01 
32 | 113,414,626 | 0,1903 0,8097 2.79E-02 1.46E-01 
33 | 117,073,162 |0,1804 0,8196 2.64E-02 1.46E-01 
34 | 120,731,698 | 0,171 0,829 2.50E-02 1.46E-01 
35 | 124,390,234 |0,1621 0,8379 2.37E-02 1.46E-01 
36 | 128,048,771 |0,1536 0,8464 2.25E-02 1.46E-01 
37 | 131,707,307 |0,1456 0,8544 2.13E-02 1.46E-01 
38 | 135,365,843 | 0,138 0,862 2.02E-02 1.46E-01 
39 /139024,38 |0,1309 0,8691 1.92E-03 1.46E-01 
40 | 142,682,916 | 0,124 0,876 1.82E-02 1.46E-01 
41 | 146,341,452 |0,1176 0,8824 1.72E-02 1.46E-01 
42 | 149,999,989 |0,1114 0,8886 1.63E-03 1.46E-01 
43 | 153,658,525 | 0,1056 0,8944 1.55E-03 1.46E-01 
44 | 157,317,061 |0,1001 0,8999 1.47E-02 1.46E-01 
45 | 160,975,598 | 0,0949 0,9051 1.39E-02 1.46E-01 
46 | 164,634,134 | 0,09 0,91 1.32E-02 1.46E-01 
47 |168292,67 |0,0853 0,9147 1.25E-02 1.46E-01 
48 | 171,951,206 | 0,0808 0,9192 1.18E-03 1.46E-01 
49 | 175,609,743 | 0,0766 0,9234 1.12E-02 1.46E-01 
50 | 179,268,279 |0,0726 0,9274 1.06E-02 1.46E-01 
51 | 182,926,815 | 0,0688 0,9312 1.01E-02 1.46E-01 
52 | 186,585,352 | 0,0653 0,9347 9.55E-03 1.46E-01 
53 | 190,243,888 | 0,0619 0,9381 9.05E-03 1.46E-01 
54 | 193,902,424 | 0,0586 0,9414 8.58E-04 1.46E-01 
55 | 197,560,961 |0,0556 0,9444 8.13E-03 1.46E-01 
56 | 201,219,497 |0,0527 0,9473 7,71E-07 1.46E-01 
57 | 204,878,033 | 0,0499 0,9501 7.31E-03 1.46E-01 
58 |208536,57 |0,0473 0,9527 6.93E-03 1.46E-01 
59 | 212,195,106 | 0,0449 0,9551 6.57E-03 1.46E-01 
60 | 215,853,642 | 0,0425 0,9575 6.22E-03 1.46E-01 
61 |219,512,178 | 0,0403 0,9597 5.90E-03 1.46E-01 
62 |223,170,715 | 0,0382 0,9618 5.59E-03 1.46E-01 


63 


226,829,251 


0,0362 


0,9638 


5.30E-03 


1.46E-01 


64 | 230,487,787 | 0,0343 0,9657 5.02E-03 1.46E-01 
65 | 234,146,324 | 0,0325 0,9675 4.76E-03 1.46E-01 
66 |237804,86 |0,0308 0,9692 4.51E-03 1.46E-01 
67 | 241,463,396 | 0,0292 0,9708 4.28E-03 1.46E-01 
68 | 245,121,933 | 0,0277 0,9723 4.06E-03 1.46E-01 
69 | 248,780,469 | 0,0263 0,9737 3.85E-04 1.46E-01 
70 | 252,439,005 | 0,0249 0,9751 3.64E-03 1.46E-01 
71 | 256,097,542 | 0,0236 0,9764 3.45E-03 1.46E-01 
72 | 259,756,078 | 0,0224 0,9776 3.27E-03 1.46E-01 
73 | 263,414,614 | 0,0212 0,9788 3.10E-04 1.46E-01 
74|267073,15 |0,0201 0,9799 2.94E-03 1.46E-01 
75 | 270,731,687 | 0,0191 0,9809 2.79E-04 1.46E-01 
76 | 274,390,223 | 0,0181 0,9819 2.64E-03 1.46E-01 
77 | 278,048,759 | 0,0171 0,9829 2.51E-03 1.46E-01 
78 | 281,707,296 | 0,0162 0,9838 2.38E-03 1.46E-01 
79 | 285,365,832 | 0,0154 0,9846 2.25E-03 1.46E-01 
80 | 289,024,368 | 0,0146 0,9854 2.13E-03 1.46E-01 
81 | 292,682,905 | 0,0138 0,9862 2.02E-03 1.46E-01 
82 | 296,341,441 | 0,0131 0,9869 1.92E-03 1.46E-01 
83 | 299,999,977 | 0,0124 0,9876 1.82E-03 1.46E-01 
84 | 303,658,514 | 0,0118 0,9882 1.72E-03 1.46E-01 
85 |307317,05 |0,0112 0,9888 1.63E-03 1.46E-01 
86 | 310,975,586 | 0,0106 0,9894 1.55E-04 1.46E-01 
87 | 314,634,122 |0,01 0,99 1.47E-03 1.46E-01 
88 | 318,292,659 | 0,0095 0,9905 1.39E-03 1.46E-01 
89 | 321,951,195 | 0,009 0,991 1.32E-03 1.46E-01 
90 | 325,609,731 | 0,0085 0,9915 1.25E-03 1.46E-01 
91 | 329,268,268 | 0,0081 0,9919 1.18E-03 1.46E-01 
92 | 332,926,804 | 0,0077 0,9923 1.12E-03 1.46E-01 
93|336585,34 |0,0073 0,9927 1.06E-03 1.46E-01 
94 | 340,243,877 | 0,0069 0,9931 1.01E-03 1.46E-01 
95 | 343,902,413 | 0,0065 0,9935 9.56E-04 1.46E-01 
96 | 347,560,949 | 0,0062 0,9938 9.06E-04 1.46E-01 


97 


351,219,486 


0,0059 


0,9941 


8.59E-04 


1.46E-01 


Após os dados mostrados será analisado as informações pelo método dos mínimos quadrados. Sabe-se 
que as falhas têm uma distribuição que segue o modelo de Weibull, desta forma, determina-se a expressão linear 
para a função cumulativa de falha. Assim, determinou-se através do Software Microsoft Excel os parâmetros 


chave para a regressão linear pelo método do rank mediano de Bernard 


4. MÉTODO DE MÍNIMOS QUADRADOS 


O método de mínimos quadrados é usualmente aplicado na análise de dados completos até a falha. A 
maior parte das distribuições aplicadas em estudos de confiabilidade podem ter suas funções facilmente 


linearizadas. Para uma distribuição de Weibull, temos as seguintes fórmulas: A confiabilidade é: 


Como verificado na técnica de verossimilhança (item 3), a distribuição que mais se a ajusta aos dados de 
tempo entre falhas é a de Weibull. Logo, aplicou-se a regressão linear por mínimos quadrados para obtenção do 


B em, visando uma posterior comparaçãoentre os valores fornecidos pelo ProConf e Excel. 


Sabendo que a confiabilidade na distribuição de Weibull é dada por: 


R(t) = PR 


Onde: = Fator de forma 
e m: Vida caracteristica 


to: Vida mínima (No nosso caso, considerada igual a “0”) 
e 1: Periodo de tempo de operação 


A função cumulativa de falhas: 


tato" 
F(t)=1-e NT 
A função densidade de probabilidade: 


“Bt-tyf o (ttf 
Ho = TS) sem 
A função taxa de falha: 


a(t) = -. (= 


Na forma linear, tem-se: 


In(—In[1 — F(t))) = Bin(t — to) — Bln(n) 


Logo, para se obter a forma linear (Y = aX + b) é necessário utilizar as seguintes 


equações: 
Y = In(=In[1 — F(to))]) 


X = In(ti—to) 


n>'P MY- MP Yi 
i=1 i=1 i=1 


a= B E I= 
nã X2- O xo? 
DD é UR DU é 
b=Pin(m) =“ cada 
n n 
cj 
n=eB 


Percebe-se, através das equações acima, que é preciso da função cumulativa pararealização dos cálculos. 
A F(t;) foi, de início, determinada pelo rank de Bernard. O B, n e TMEF foram adquiridos pela estimativa 
mostrada nas tabelas a seguir. A confiabilidadeintrínseca tO foi assumida com o valor 0, 

Para obtermos os parâmetros das equações acima, usaremos a regressão linear por Mínimos Quadrados. 
A partir da equação da confiabilidade, ao aplicarmos o logaritmo natural nos dois lados da equação na forma 


linear, e deduzindo chegaremos a: 


528.53 | 0.189266 6.270107 | -1.56153 2438378 | 39.31424 | -9.79096 
552.53 | 0.217514 6.314501 | -1.40536 1.975029 | 39.87293 | -8.87413 
576.18 | 0.245763 6.356428 | -1.26568 1.601939 | 40.40418| -8.04519 
597.36 | 0.274011 6.392528 | -1.13874 1.296739 | 40.86441 | -7.27946 
624.61 | 0.30226 6.437129 | -1.02191 1.04429| 41.43663 | -6.57814 
671.61 | 0.330508 6.509683 -0.9132 0.833941 | 42.37597 | -5.94467 
671.79 | 0.358757 6.509941 | -0.81115 0.657964 | 42.37933 | -5.28054 
671.94 | 0.387006 6.510166 | -0.71458 0.510619| 42.38227| -4.65201 
672.29 | 0.415254 6.510693 | -0.62254 0.38756| 42.38913| -4.05319 
695.71 | 0.443503 6.544928 | -0.53428 0.285451| 42.83608| -3.4968 
695.86 | 0.471751 6.545146| -0.44912 0.20171| 42.83894| -2.93957 
GuIBaSo! 0.5 6.576921 | -0.36651 0.134332 | 43.25589| -2.41053 
718.98 | 0.528249 6.57784 | -0.28595 0.081765 | 43.26798| -1.88091 
719.64 | 0.556497 6.578747 | -0.20696 0.042833 | 43.27991| -1.36155 
719.88 | 0.584746 6.57908 | -0.12912 0.016673 43.2843 | -0.84952 
720.41 | 0.612994 6.579816 | -0.05201 0.002705 | 43.29398 | -0.34224 
746.57 | 0.641243 6.615493 0.0248 0.000615 | 43.76475 | 0.164063 
768.18 | 0.669492 6.644023 | 0.101765 0.010356| 44.14304 | 0.676127 
768.49 | 0.69774 6.644424 | 0.179374 0.032175 | 44.14837/1.191837 
768.84 | 0.725989 6.644884 | 0.258191 0.066663 | 44.15449 | 1.715649 
791.30 | 0.754237 6.673671 0.33889 0.114846| 44.53788| 2.26164 
792.07 | 0.782486 6.674655 | 0.422317 0.178351| 44.55102 |2.818818 
814.99 | 0.810734 6.70317 | 0.509587 0.259679 | 44.93249|3.415851 
816.10 | 0.838983 6.704541 | 0.602262 0.36272 | 44.95087 | 4.037892 
816.31 | 0.867232 6.704794 | 0.702676 0.493754| 44.95426| 4.7113 
839.91 | 0.89548 6.733292 | 0.814647 0.66365 | 45.33722 | 5.485258 
839.91 | 0.923729 6.733299 | 0.945251 0.8935 | 45.33731 | 6.364659 
840.18 | 0.951977 6.733619| 1.110568 1233362 | 45.34163] 7.478145 
840.49 | 0.980226 6.733984 | 1.366955 1.868567| 45.34654 | 9.205054 
SOMA 227.8633 | -19.4526 60.1767 | 1484.726 | -118.946 


Tabela 5 - Valores estimados fundamentais para determinação do B, 7 e TMEF. 


3.4147 


b -19.224 
B 3.4147 
n 278.9507 


Tabela 6 — Parâmetros de Weibull 


Podemos comparar o método dos mínimos quadrados como o método daverossimilhança pelo valor 


de TMEF, através dos novos parâmetros encontrados. Logo: 


TMEF =0 + 73284,967x T(1 + 1/7,0775) 


TMEF = 73284,967 x T(1,14) 


TMEF = 73284,967x (1,14) 


TMEF = 685,95 horasTabela 4Cálculo da confiabilidade, função cumulativa, densidade 
de falha e taxa de falha pra distribuição Weibull. 


Densidade de 


índices |t Expoente Weibull | Confiabilidade | Função cumulativa (F) | Taxa de falha (À) | falha (f) 
1 431.59375 -4.438625478 0.0118123 0.988187837 0.035118 0.000415 
2 |455.6409722 -5.341418051 0.0047895 0.995210925 0.04003 0.000192 
3 |503.2777778 -7.501020234 0.000552 0.99944748 0.050894 2.81E-05 
4 |504.1923611 -7.547669096 0.000527338 0.999472662 0.051117 2.7E-05 
5 |504.2083333 -7.548485587 0.000526907 0.999473093 0.051121 2.69E-05 
6 |528.2451389 -8.849603456 0.000143439 0.999856561 0.057206 8.21E-06 
7 |528.5340278 -8.866140521 0.000141086 0.999858914 0.057281 8.08E-06 
8 |552.5263889 -10.31740008 3.30529E-05 0.999966947 0.063763 2.11E-06 
9 |576.1847222 -11.905498 6.75317E-06 0.999993247 0.070557 4.76E-07 
10 |597.3645833 -13.46735528 1.41645E-06 0.999998584 0.076983 1.09E-07 
11 | 624.6111111 -15.68288143 1.54529E-07 0.999999845 0.085737 1.32E-08 
12 | 671.6131944 -20.09191799 1.88014E-09 0.999999998 0.102154 1.92E-10 
13 | 671.7868056 -20.10965857 1.84708E-09 0.999999998 0.102218 1.89E-10 
14 | 671.9381944 -20.12513739 1.81871E-09 0.999999998 0.102273 1.86E-10 
15 | 672.2923611 -20.16138223 1.75397E-09 0.999999998 0.102403 1.8E-10 
16 695.70625 -22.66152515 1.43954E-10 1 0.111228 1.6E-11 
17 | 695.8583333 -22.6784456 1.41539E-10 1 0.111287 1.58E-11 
18 | 718.3243056 -25.27756404 1.05219E-11 1 0.120162 1.26E-12 
19 | 718.9847222 -25.35 700917 9.71831E-12 1 0.120429 1.17E-12 
20 | 719.6368056 -25.43562491 8.98356E-12 1 0.120693 1.08E-12 
21/719.8770833 -25.46463642 8./2668E-12 1 0.12079 1.05E-12 
22 | 720.4069444 -25.52869534 8.18518E-12 1 0.121005 9.9E-13 
23 | 746.5729167 -28.83614129 2.99656E-13 1 0.131892 3.95E-14 
24 | 768.17/91667 -31.78676131 1.56742E-14 1 0.141298 2.21E-15 
25 768.4875 -31.83034934 1.50057E-14 1 0.141435 2.12E-15 
26 | 768.8409722 -31.88037048 1.42735E-14 1 0.141592 2.02E-15 
27 | 791.2951389 -35.1733607 5.30156E-16 1 0.151785 8.05E-17 
28 | 792.0743056 -35.29176693 4.70956E-16 1 0.152146 7.17E-17 
29 | 814.9854167 -38.90099841 1.275E-17 1 0.162991 2.08E-18 
30 | 816.1034722 -39.08353336 1.06228E-17 1 0.163531 1.74E-18 
31 | 816.3097222 -39.11727201 1.02703E-17 1 0.163631 1.68E-18 
32 | 839.9076389 -43.1152399 1.8849E-19 1 0.175288 3.3E-20 


33 | 839.9131944 -43.11621373 1.88307E-19 1 0.175291 3.3E-20 


34 | 840.1826389 -43.16346312 1.79616E-19 1 0.175426 3.15E-20 


35 | 840.4888889 -43.21721112 1.70217E-19 1 0 


Tabela 7 - Cálculo da confiabilidade, função cumulativa, densidade de falha e taxa de falha pra 


distribuição Weibull. 


Confiabilidade 
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Conforme previsto, o ajuste das curvas à distribuição Weibull foi capaz de ajustaras geometrias de tal 
forma a favorecer a compreensão do comportamento das funções. A taxa de falha, evidentemente não linear 
(como supõe a distribuição exponencial), tende aaumentar, característica típica de sistemas que, assim como o 
equipamento hipotético aqui abordado, tendem a aumentar sua probabilidade de falhas, devido a fatores diversos, 
como corrosão, fadiga e desgaste. Isso fica ainda mais evidente quando se observa o valordo fator de forma p = 
3,73 >1, cuja consequência se reflete em características de falha por desgaste, com mortalidade infantil 
inexpressiva. Admitindo um crescimento aproximadamente linear para A(t), portanto, o equipamento em questão 
situa-se entre o grupo dos 5% dos componentes que apresentam essa característica. 

Verifica-se que, apesar da aparente complexidade das equações de linearização, os resultados se aproximam de 
distribuições mais plausíveis de dados. Além disso, os parâmetros obtidos podem ser utilizados posteriormente, 


na análise de custos. 


4.2 - REPARO 
Como verificado na técnica de verossimilhança (item 3), a distribuição que maisse a ajusta aos 


dados de tempos de reparo é a Lognormal. Para este tipo de distribuição,a função densidade de 


probabilidade é dada por: 


mt) = pasa e 


A função distribuição cumulativa é expressa por: 


Ntij= d [Ei Ê did 


O tempo médio para reparo é dado por: 


2 


TMPR = Be 7 


Logo, a forma linear é expressa por: 
InCti) = ap *[MCt)] + In(B) 
Onde: 
Y= In(t;) 


on Bisa Xi Vi — Bisa Mi Dia 
nãE 1 x2 “Er sã) 


X = q" [MCt)] 


ny a 
Es In(B) = Di=1 t no Rx t 
n n 


Sendo a o parâmetro de forma ou dispersão (fornece a variabilidade das equipes de 
manutenção), B o parâmetro de escala (fornece o desempenho das equipes de manutenção), q a 


função normal padronizada e "ta inversa da função normal padronizada. 


Percebe-se, através das equações acima, que é preciso da função mantenabilidade para 
realização dos cálculos. A M(ti) foi, de início, determinada pelo rank de Bernard. Oq, B e TMPR 


foram adquiridos pela estimativa mostrada nas tabelas 10 e 11. 


EE 
1 2.1 0.028571429 | 0.741937 | -1.902216496 | 0.550471 | 3.618428 | -1.41133 
Rossas 0.057142857 | 0.96825| -1.579219519 | 0.937509 | 2.493934 | -1.52908 
3 /3.433333 0.085714286 | 1.233532 | -1.367627923 | 1.5216] 1.870406 | -1.68701 
ao us 0.114285714 | 1.262242 -1.20404696 | 1.593254 | 1.449729 | -1.5198 
5 |4.166667 0.142857143 | 1.427116| -1.067570524 | 2.036661 | 1.139707 | -1.52355 
6 |4.833333 0.171428571 | 1.575536| -0.948535041 | 2.482315 | 0.899719 | -1.49445 
pets 0.2/1.664557| -0.841621234 | 2.770751 | 0.708326 | -1.40093 
8 5.65 0.228571429 | 1.731656| -0.743559757 | 2.998631 | 0.552881 | -1.28759 
9 bi! 0.257142857| 1.7492| -0.652178986| 3.0597|0.425337 | -1.14079 
10 6.4 0.285714286 | 1.856298 | -0.565948822 | 3.445842 | 0.320298 | -1.05057 
11/6.533333 0.314285714 | 1.876917 | -0.483738553 | 3.522818 | 0.234003 | -0.90794 
12 6.95 0.342857143 | 1.938742 | -0.404677905 | 3.758719 | 0.163764 | -0.78457 
13 | 6.966667 0.371428571 | 1.941137 | -0.328072108 | 3.768012 | 0.107631 | -0.63683 
14 725 0.4/1.981001| -0.253347103 | 3.924367 | 0.064185 | -0.50188 
15/7.316667 0.428571429 | 1.990155| - -0.18001237 | 3.960716 | 0.032404 | -0.35825 
16 76 0.457142857 | 2.028148 | -0.107634392 |4.113385 | 0.011585| -0.2183 
17 | 8.283333 0.485714286 | 2.114245 | -0.035816632 | 4.470034 | 0.001283 | -0.07573 
18/8.316667 0.514285714 |2.118262 | 0.035816632 | 4.487032 | 0.001283 | 0.075869 
19 8.6 0.542857143 |2.151762 |  0.107634392 | 4.630081 | 0.011585 | 0.231604 
20 | 8.783333 0.571428571 | 2.172856 0.18001237 | 4.721303 | 0.032404 | 0.391141 
2 | 8783335 0.6 |2.172856| | 0.253347103 | 4.721303 | 0.064185 | 0.550487 
Pupo 8.95 0.628571429 |2.191654|  0.328072108 | 4.803345 | 0.107631 | 0.71902 
23 |9.216667 0.657142857 |2.221013 |  0.404677905 | 4.932901 | 0.163764 | 0.898795 
24 |9.733333 0.685714286 |2.275556|  0.483738553 | 5.178157 | 0.234003 | 1.100774 


us as 0.714285714 |2.277267 | 0.565948822 | 5.185946 | 0.320298 | 1.288817 
26 9.85 0.742857143 |2.287471 |  0.652178986 | 5.232526 | 0.425337 | 1.491841 
Bi 0/5833 0.771428571 | 2.35928| 0.743559757 | 5.566204 | 0.552881 | 1.754266 
28 | 11.58333 0.8 |2.449567 |  0.841621234 | 6.00038 | 0.708326 | 2.061608 
29 11.65 0.828571429 |2.455306| 0.948535041 | 6.028528 | 0.899719 | 2.328944 
30 | 12.33333 0.857142857 |2.512306| 1.067570524| 6.31168| 1.139707 | 2.682063 
31/13.16667 0.885714286 | 2.577688 1.20404696 | 6.644477 | 1.449729 | 3.103658 
a” ISOS 0.914285714 | 2.63548| 1.367627923 | 6.945752 | 1.870406 | 3.604355 
33 | 17.31667 0.942857143 |2.851669| 1.579219519| 8.132019 | 2.493934 | 4.503412 
34 18.55 0.971428571 | 2.92047| 1.902216496 | 8.529144 | 3.618428 | 5.555366 
aS O OB3aS 1/2.994899 8.969418 

SOMA 71.70603 | 1.15463E-14| 155.935 | 28.18724 | 14.81343 


Tabela 8 - Valores estimados fundamentais para determinação do a, B e TMPR. 


Tabela 5 — Parâmetros calculados 


a [052554] 

bo |2o4874| 

o (0.52554| 

Po Jpisso 

Ep] 

ju [890701] 

[media — [2.04874] 
24| — 0.006448291 24| 0.016425702 24 0.016532307 
267] 0.019892215| 267| 0.034835173| 2.67 0.035542186 
395| — 0.060427011| 3.95] 0.066404485| 3.95 0.070675174 
425| 0067252738) 4.25] 0070126752) 4.25 0.075183016 
4.28 0.11843574] 4.28| 0090534448] 4.28 0.102697502 
443| 0183946592]  443| 0104740379) 4.43 0.128349907 


4.5 0.232378444 4.5 | 0.110017659 4.5 0.143322784 
4.6 0.273134393 4.6| 0.112025627 4.6 0.154121512 
4.98 0.284346627 4.98 | 0.112254483 4.98 0.15685594 
5.03 0.357112022 5.03 | 0.110947826 5.03 0.17257723 
Sodl 0.371851227 5.1| 0.110171518 5.1 0.175390804 
5.15 0.417101614 5.15 | 0.106885422 Sodls) 0.183368876 
5.35 0.418881297 5.35 | 0.106730378 5.35 0.183663642 
5.37 0.448717906 5.37 | 0.103865419 5.37 0.188407024 
5.37 0.45561627 5.37 | 0.103134715 5.37 0.189452237 
5.43 0.484369625 5.43 | 0.099832031 5.43 0.193611618 
5.68 0.549594771 5.68|  0.090957514 5.68 0.201945954 
6.02 0.552618367 6.02 | 0.090504067 6.02 0.202297233 
6.87 0.577704724 6.87 | 0.086611286 6.87 0.205096506 
7.15 0.593346769 7.15 | 0.084071099 7.15 0.206739041 
7.25 0.593346769 7.25 | 0.084071099 225) 0.206739041 
8.32 0.607162433 8.32 | 0.081759213 8.32 0.20812473 
8.47 0.628467583 8.47| 0.078074713 8.47 0.210142398 
9.05 0.666977629 9.05 | 0.071073675 9.05 0.213420122 
9.1 0.668160044 9.1, 0.070852188 9.1 0.213513132 
9.12 0.675177281 9.12 | 0.069530122 ORIZ! 0.214055601 
9.15 0.722704309 9.15 0.06025267 S)ál5) 0.217286717 
9.23 0.777176855 9.23| 0.049008188 9.23 0.219942087 
9.77 0.780420286 9.77 | 0.048320707 9.77 0.22005998 
10.15 0.811131304 10.15 | 0.041723868 10.15 0.22091468 
10.28 0.84290866 10.28 |  0.034751112 10.28 0.221215964 
10.47 0.86788525 10.47 | 0.029186179 10.47 0.22091537 
10.92 0.936721358 10.92 | 0.013648357 10.92 0.215686628 
2572 0.951415322 12.52 | 0.010342233 IRS) 0.21287026 
13.43 0.964097981 13.43 0.0075147 13.43 0.20931136 


Tabela 11 - Cálculo da Manutenibilidade, função densidade m(t), densidade de falha etaxa de falha para 


distribuição Lognormal com os q, B e TMPR determinados anteriormente. 


Tabela 12 - Cálculo da Manutenibilidade, função densidade m(t), densidade de falha e taxa de falha para distribuição Lognormal com os a, 8 e 


TMPR determinados anteriormente. 


Taxa de 
Manutebalidade acumulativa | Densidade Falha 
0.006448291 0.016425702 | 0.016532307 
0.019892215 0.034835173 | 0.035542186 
0.060427011 0.066404485 | 0.070675174 
0.067252738 0.070126752 | 0.075183016 
0.11843574 0.090534448 | 0.102697502 


0.183946592 


0.104740379 


0.128349907 


0.232378444 0.110017659 | 0.143322784 
0.273134393 0.112025627 | 0.154121512 
0.284346627 0.112254483 | 0.15685594 
0.357112022 0.110947826 | 0.17257723 


0.371851227 


0.110171518 


0.175390804 


0.417101614 0.106885422 | 0.183368876 
0.418881297 0.106730378 | 0.183663642 
0.448717906 0.103865419 | 0.188407024 

0.45561627 0.103134715 | 0.189452237 
0.484369625 0.099832031 | 0.193611618 


0.549594771 


0.090957514 


0.201945954 


0.552618367 0.090504067 | 0.202297233 
0.577704724 0.086611286 | 0.205096506 
0.593346769 0.084071099 | 0.206739041 
0.593346769 0.084071099 | 0.206739041 
0.607162433 0.081759213 | 0.20812473 
0.628467583 0.078074713 | 0.210142398 
0.666977629 0.071073675 | 0.213420122 
0.668160044 0.070852188 | 0.213513132 
0.675177281 0.069530122 | 0.214055601 
0.722704309 0.06025267 | 0.217286717 


0.777176855 0.049008188 | 0.219942087 
0.780420286 0.048320707 | 0.22005998 
0.811131304 0.041723868 | 0.22091468 

0.84290866 0.034751112 | 0.221215964 

0.86788525 0.029186179 | 0.22091537 
0.936721358 0.013648357 | 0.215686628 
0.951415322 0.010342233 | 0.21287026 
0.964097981 0.0075147 | 0.20931136 


Manutebalidade acumulativa 


Manutebaidade acumulativa 


Figura 33-Gráfico Função Manutenabilidade Cumulativa - 


Densidade Figura 34 - Gráfico Função 


Densidade de Probabilidade 


O resultado obtido pelo cálculo dos parâmetros da distribuição lognormal aos tempos dereparo atingiu 
resultados semelhantes àqueles verificados para a linearização das equações de Weibull para dados de falha, 
seção 4.1. Verifica-se um comportamento crescente na taxa de reparo, o que é esperado para um componente 
cuja característica defalha é crescente, conforme verificado anteriormente naquela seção. 

Uma vez que a máquina tende a falhar com maior frequência, as paradas para reparo serão naturalmente 
mais frequentes. O método dos mínimos quadrados é, em geral, menos eficaz que o método da verossimilhança, 


desenvolvido na seção 3. Apesar de ambos objetivarem estimar parâmetros para adequar o conjunto de dados 


de falha/reparo às distribuições padrão, os cálculos mais complexos - e comumente executados por programas 
específicos — do segundo obtém resultados que conferem melhores características, como: não tendenciosidade, 
consistência, eficiência e suficiência. 
5- Testes de Aderência 
Os resultados mostrados na seção 3 permitem uma avaliação rápida e qualitativada 
possibilidade de ajuste dos dados com os modelos matemáticos ideais. A validação ou não 
desses modelos como representativos da amostra depende do resultado dos Testes de 
Hipótese, a partir da seguinte premissa: a hipótese de que o modelo é representativo da 
amostra é denominada hipótese nula H$, ao passo que a afirmação oposta é denominada 
H1. Os testes são de natureza quantitativa, e, quando não analisados corretamente, podem 


ocasionar erros, conforme ilustrado na tabela 13. 


Teste de ipótese 
Ho verdadeiro H, verdadeiro 
Aceitar Ho Decisão correta Erro Tipo IH 
Aceitar Hi Erro Tipo I Decisão correta 
1 


Os testes realizados nesse estudo utilizaram tanto a ferramenta do ProConf para testes deaderência, quanto 


o cálculo em planilha eletrônica, para fins de comparação. 


5.1.5.1- Teste do Qui-quadrado (y?) 


O teste do qui-quadrado consiste na subdivisão dos dados em classes, e em seguida, na 
estimativa do grau de aproximação entre a frequência de dados em um determinado intervalo e 


a frequência esperada — estatisticamente — para aquele intervalo. Matematicamente: 


k 
(Oi — Ei)? 
Ei 


É, 
| 
M 


k = número de classes (pela regra de Sturges, k = 1 + 


3,3log10n; 


Oi = número de falhas (ou reparos) observados na 1- 


éssima classe; 


Ei = número de falhas (ou reparos) esperados na 1- 


éssima classe; 


n = número total em risco (tamanho da amostra); 


pi = probabilidade da falha ocorrer na i-éssima classe, se HO é verdadeira; 


A esperança é dada por: 


E =nPi 
Tal que: 


Pi=F(ai) — F(ai1) 


O processamento dos dados no ProConf fornece a estimativa do 


parâmetro x2. Esse encontra-se tabelado para seus valores críticos (acima dos 
quais a hipótese nula é rejeitada) níveis diferentes de significância e graus de 


liberdade diferentes. 


5.2 — Teste Kolmogorov-Smirnov 


O teste possui a mesma finalidade do anterior, e, para uma dada amostra de dados 
agrupados, consiste no cálculo — e posterior comparação com os valores críticos tabelados 


— do parâmetro D, expresso por: 


Onde, 


k = número de classes; 
Ox = número de falhas (ou reparos) observados na 1- 
éssima classe; 


Ex = número de falhas (ou reparos) esperados na 1-éssima 


classe; 


N = número total em risco (tamanho da amostra); 


pi = probabilidade da falha ocorrer na i-éssima classe, se H, é verdadeira; 


5.3 - Resultado dos testes de aderência 


Resultado a partir do Proconf: 


O processamento dos dados forneceu como resultado dos dois testes — realizados 
simultaneamente pelo software. A figura 23 é referente aos dados de falha do 


equipamento hipotético. 


A esperança é dada por: 


Eij=nPi 


Tal que: 


Pi; = F(ai) — F(az1) 


O processamento dos dados no ProConf fornece a estimativa do parâmetro 2. Esse 
encontra-se tabelado para seus valores críticos (acima dos quais a hipótese nula é rejeitada) 


níveis diferentes de significância e graus de liberdade diferentes. 


5.4 — Teste Kolmogorov-Smirnov — (K-S) 


O teste possui a mesma finalidade do anterior, e, para uma dada amostra de dados 
agrupados, consiste no cálculo — e posterior comparação com os valores críticos tabelados 


— do parâmetro D, expresso por: 


Ox — Ek 
Ek 


D = max | | < Dor 


k = número de classes; 


Ox = número de falhas (ou reparos) observados na i-éssima classe; 
Ex = número de falhas (ou reparos) esperados na 1-éssima classe; 
N = número total em risco (tamanho da amostra); 


pi = probabilidade da falha ocorrer na i-éssima classe, se HO é verdadeira; 


5.5 - RESULTADO DOS TESTES DE ADERÊNCIA POR PRONCOF 


O processamento dos dados forneceu como resultado dos dois testes — realizados 


simultaneamente pelo software. A figura 23 é referente aos dados de falha do equipamento 


hipotético. 


“S Ajuste do Modelo de Distribuição EEIS 
Figura 35 - Teste de 
Modelos | Ajuste/E statísticas | Funções de Confiabilidade | Gráficos Testes de Aderência | Aderência por Proncof 


Teste Analítico do 72: 


or : Ee : Fregiiência Frequência Valor do 
Limite Inferior | Limite Superior Dsaveda qui-quadrado 


41114,5 0,01 


Z2= 5.28 com 3 graus de liberdade Nível de Significância = 0,1523 
Teste Analítico de Kolmogorov-Smirnov (K-5): 
DN = 0,1149 Nível de Significância = 0,1881 


A hipótese de que a população segue o modelo Weibull 
não pode ser rejeitada. 


[E &) o de 
1 


Observa-se que o tamanho dos intervalos não obedece à regra estabelecida para 
o teste qui-quadrado. A distribuição assumida para a hipótese nula foi a de Weibull. Os 


valores de x2 são somados para cada intervalo, totalizando 7,82. De acordo com a tabela 


disponível na literatura, o valor crítico desse parâmetro para significância a = 0,10 e 3 


graus de liberdade é 4,87. Portanto, y? > xcr 2, logo, a hipótese HO pode ser rejeitada 


de acordo com o teste qui-quadrado. 


O valor calculado para D (DN) é de 0,13. De acordo com a literatura, o teste K- 
S fornece resultados mais precisos que o anterior, o que o torna ferramenta essencial 
para se tomar conclusões acerca do modelo estatístico condizente com a amostra. Na 
literatura, encontra-se que Dcrt. = 0,202 para um tamanho de amostra N=35 e 
significância 
a = 0,10. Logo, D < Dcrt, e a amostra apresenta-se verossímil ao modelo 
Weibull. Essepadrão de falhas é comum para boa parte dos componentes 


mecânicos. 


Resultados similares foram obtidos para os tempos de reparo, conforme a 
figura 36 


“Ajuste do Modelo de Distribuição [osPa Figura 36 - teste 


E E: = E: E PEER de aderencia 
Modelos | Ajuste/Estatísticas | Funções de Confiabilidade | Gráficos Testes de Aderência | 


Teste Analítico do X2: 


e : a : Frequência Frequência Valor do 
Limite Inferior | Limite Superior DEsSveda qui-quadrado 


411145 

56107,5 
561075 

711005 
711005 


Z2= 5.69 com 3 graus de liberdade Nível de Significância = 0,1279 
Teste Analítico de Kolmogorov-Smirnov (K-5): 
DN = 0,1842 Nível de Significância = 0,0041 


A hipótese de que a população segue o modelo Lognormal 
é rejeitada. 


[E] &) | e, 


“N Ajuste do Modelo de Distribuição ooo 


Modelos | Ajuste/E statísticas | Funções de Confiabilidade | Gráficos Testes de Aderência | 


Teste Analítico do X2: 


as E nã : Freguência Frequência Valor do 
Limite Inferior | Limite Superior Obsersada qui-quadrado 


411145 5 016 
561075 E 


56107,5 
711005 
711005 


%2= 5.69 com 3 graus de liberdade Nível de Significância = 0,1279 


Teste Analítico de Kolmogorov-Smirnov (K-5): 
DN = 0,1842 Nível de Significância = 0,0041 


A hipótese de que a população segue o modelo Lognormal 
é rejeitada. 


[E] & aj de, 


Apesar dos dois resultados, verifica-se que o primeiro caso apresentou 
parâmetrosmais distantes dos valores críticos — para essas condições de significância, 
graus de liberdade e tamanho de amostra — ou seja: Xiognormal < Xgamma < Xcrt € 
Dgamma < Diognormal < Dart. 

Assim, a distribuição Lognormal é mais representativa da amostra que a 
distribuição gama. Esse resultado também era esperado, uma vez que ocomportamento 
dos tempos de reparo apresenta variabilidade em função de fatores humanos. Nesse 
sentido, os tempos de reparo podem variam muito de acordo com o nívelde treinamento, 
facilidade de diagnose, curva de aprendizado típica do ser humano, edemais fatores 


humanos e organizacionais. 


5.6 — Metodos pelo Excel até a falha 


Os testes de aderência anteriores foram executados na planilha eletrônica, para fins de 


comprovação. Considerando uma amostra com 35 tempos de falha, obtemos um k=6 erealizamos o teste 


de aderência Qui-Quadrado para os dados de falha: 


1 125,95 5 

13 0,130387682 4,56 15,595989 
2 292,59 8 
3 378,37 6 

8) 0,198784528 6,96 0,5996408 
4 464,15 3 
5 549,93 10 9 0,23093836 8,08 0,1040695 
6 635,71 3 S 0,261801817 9,16 0,0029018 


Tabela 18 - Teste de aderência Qui-Quadrado para os dados de falha 


Com uma amostra de 35 tempo de falha, obtemos um k=6 e realizamoso teste de 


aderência Kolmogorov-Smirnov: 


Tabela 19 - Teste de aderência Kolmogorov-Smirnov para os dados de falha 


5.7 — Reparo dos métodos 


Tabela 20 - Teste de aderência Qui-Quadrado para os dados de falha 


Considerando uma amostra de 35 tempo de falha, obtemos um k=6 e realizamos 


o teste de aderência Kolmogorov-Smirnov: 


Somatório: 0,50888 


Tabela 21 - Teste de aderência Kolmogorov-Smirnov para os dados de falha 


6- custos de manutenção preventiva X manutenção corretiva 


Como a distribuição do tempo entre falhas segue o modelo de Weibull, a equaçãoque 
compara o custo com manutenção preventiva (Cmp) e com manutenção corretiva (Cmc) é 


expressa por: 


(Cc) É expressa por: 


(mp  U—to (O 
Cia E = Ea es 


E o desvio padrão é determinado por: 


o =" +28) =D + pj? 


Sendo o TMEF, Kc o custo de uma intervenção de manutenção corretiva e Kp o custo 
de uma intervenção de manutenção preventiva. Para valores de taxas de custos Kp/Kc iguais 


a 0,6; 0,4; 0,2 e 0,1 obteve-se a relação Cmp/Cmc como mostra a tabela 18. 


| | 


0,720 
,119 
0,737 


0,773 
,069 0,813 
,094 0,862 


0,877 
0,936 


Tabela 22 - Relações Cmp/Cmc para valores de Kp/Kc 


Além disso, plotou-se um gráfico comparativo da relação entre o custo horário da 
manutenção preventiva e o custo horário da manutenção corretiva, mostrado a seguir. 


10,000 
9,000 


8,000 
Figura 37 - 
Relações 7,000 
Cmp/Cmc, para 


diversos valores 6,000 =E-Kp/Kc=0,6 
alia 5,000 —4—Kp/Kc=0,4 
4.000 4 Kp/Kc=0,2 

—d&— Kp/Kc=0,1 


3,000 


2,000 


1,000 


0,000 


2 92 183 275 367 457 548 640 732 822 913 1005 


A relação de custo foi plotada versus os tempos entre falhas disponíveis. Pela análise do 
gráfico, percebe-se que para Kp/Kc iguais a 0,1 e 0,2 a manutenção preventivaé mais econômica 
até um tempo de, aproximadamente, 500 h. Já para Kp/Kc iguais a 0,4e 0,6 a manutenção 


corretiva é sempre a mais econômica. 


7- PERÍODOS IDEAIS PARA REALIZAR A MANUTENÇÃO PREVENTIVA 


A determinação do período “T”, para uma política de manutenção preventiva 
sistemática, cujo objetivo seja minimizar os gastos da manutenção do equipamento ou 
componente é realizada através de duas etapas. Primeiramente, a seguinte condição deveser 


satisfeita: 


Por fim, determina-se o período padrão ótimo “T” através de: 


1/4/9 
T=to+n ARE, 


KB-1 


O tempo ótimo pode ser verificado na tabela 19, para valores de kp/kc<0,391. 


Períodos Ótimos — Por tempo 


|  Tempo(r) | 


oa 167,89 


TABELA 24: Intervalos ideais para manutenção preventiva sistemática 


7.1 — Modelo relacionado com a idade 
Nesse caso, os intervalos de tempo são definidos a partir da idade do 
componente,e não permanecem constantes, ao contrario do método da seção 7.1. 
Em um intervalo decomprimento tp, entre substituições preventivas, existem duas 


possibilidades: 


a) O componente não apresenta falhas aleatórias ao longo do intervalo; 
b) O componente apresenta falhas ao longo do intervalo. 


A probabilidade de ocorrência de cada um desses eventos está ligada às 


funções de falha, como F(t) e f(t). O custo total de um intervalo deve abranger, 


probabilisticamente, as duas possibilidades, e nesse sentido, é possível mostrar 


que: 


C(tp) = R 
Jo” RCDdt 


Para analisar o comportamento do custo total, foi utilizado o software MATLABTM, e 
os gráficos foram plotados a partir do seguinte script, no qual o valor docusto de manutenção 
preventiva foi arbitrado em R$ 1000,00 (o valor não altera qualitativamente a forma da curva, 
nem o ponto de mínimo, que é o que se busca), e o fator multiplicador de Cf foi modificado 
para 0,1, 0,2, 0,4 e 0,6. As curvas da figura 35 foram obtidas para os respectivos valores de 


razão. 


%% Cálculo do custo e tempo mínimos por idade 


Cp=1000;%custo de manuteção preventiva 
Cf1l=Cp*(1/0.1); 


Cf2=Cp*(1/0.2); 
Cf3=Cp* (1/0.4); 


Cf4-Cp* (1/0.6);%custo de manutenção corretiva 


beta=2.264;%fator de forma 


neta=514.82;%fator de escala 
mintp=10;&menor intervalo considerado 
maxtp=1000;%maior intervalo considerado 


tp=mintp:1l:maxtp;%valores de intervalo possíveis na análise 
R=exp (- (tp. /neta) .“beta); “confiabilidade no tempo tp intervalar 
CtpNUMI=Cp. *R+C£1.*(1-R); 


CtpNUM2=Cp. *R+Cf2.*(1-R); 
CtpNUM3=Cp. *R+Cf3.*(1-R); 
CtpNUM4=Cp. *R+Cf4.*(1-R); 


Fun=€ (x) exp(-(x./neta) .“beta); 
Denoml=integral (Fun, 0,mintp); 
for i=mintp+1l:1l:maxtp 


Denom2=integral (Fun,0,1); 
Denomi=[Denom1, Denom2]; 


end 


Custo1l=CtpNUM1./Denoml; 
Custo2=CtpNUM2. /Denoml; 
Custo3=CtpNUM3. /Denoml; 
Custo4=CtpNUM4. /Denoml; 
figure (1) 

subplot (2,2,1) 


plot (tp, Custol, 'linewidth',3);grid; 


xlabel ("Intervalo tp entre substituições (h)'); 


ylabel ('Custo total esperado por intervalo(RS$S)')J;ylim([0 110]); 


subplot (2,2,2) 
plot (tp, Custo2, 'linewidth',3);grid; 
xlabel ("Intervalo tp entre substituições (h)'); 


ylabel ('Custo total esperado por intervalo(R$S)');ylim([0 110]); 


subplot (2,2,3) 
plot (tp, Custo3, 'linewidth',3);grid; 
xlabel ("Intervalo tp entre substituições (h)'); 


ylabel ('Custo total esperado por intervalo(RS$S)')J;ylim([0 110]); 


subplot (2,2,4) 
plot (tp, Custo4, 'linewidth',3);grid; 
xlabel ("Intervalo tp entre substituições (h)'); 


ylabel ('Custo total esperado por intervalo(RS$S)')J;ylim([0 110]); 


[Custol, tp]=min (Custol); 
intervalol=tp+mintp; 
fprintf('Intervalo ideal 0,1:'); 
disp (intervalol); 
[Custo2, tp]=min (Custo2); 
intervalo2=tp+mintp; 
fprintf('Intervalo ideal 0,2:'); 
disp (intervalo2); 
[Custo3, tp]=min (Custo3); 
intervalo3=tp+mintp; 
fprintf('Intervalo ideal 0,4:'); 
disp (intervalo3); 
[Custo4, tp]=min (Custo4) ; 
intervalo4=tp+mintp; 
fprintf('Intervalo ideal 0,6:'); 
disp (intervalo4); 


Figura 38 - Curvas 
de custo por idade 
para diferentes 
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Os tempos onde ocorre o menor custo total para cada curva estão mostrados na tabela 25 


TABELA 25: Intervalos ideais para manutenção 
preventiva sistemática (por idade) 


8 — Método Monte Carlo 


O equipamento hipotético foi submetido a uma simulação de duração 
total de 40000 horas, considerando o intervalo ideal para manutenção 
preventiva sistemática de 228 h (Tempo ótimo obtido para Kp/Kc=0,2). 

Para cálculo do TEF foram usados os parâmetros B, n obtidos 
anteriormente (Distribuição Weibull) e para cálculo do TPR foram usados os 
parâmetros q, B log obtidos anteriormente (Distribuição Lognormal). 

O script do Matlab mostrado a seguir seguiu a seguinte linha de 
raciocínio: Como tempo ótimo estabelecido, calcula-se inicialmente o TEF. Se 
ele der menor que o tempoótimo, é calculado o TPR. Ou seja, os TEF e TPR 
são calculados até atingirem o valor de 228 horas. Quando esse valor é 
alcançado, existe um Tp (Tempo de manutenção preventiva), estabelecido 
como 0,8*TPR. Depois de contabilizado o Tp, um novo intervalo de 228 horas 
é estabelecido, para cálculo de novos TEF e TPR. Esses procedimentos são 
efetuados até ser atingido o tempo de duração da simulação adotado. No caso, 


40000 horas. 


clear all 
close all 
clc 


TMPR=10.28;%Tempo Médio para reparo 


Ts=40000; % tempo de Simulação 


Tot=228; % Intervalo ótimo 
To=0; %Tempo inicial do intervalo 
rel=0; % Tempo corrido 


cont=0.00000001;  SContador 


Tfalha=0; &Tempo de falha em To 


Treparo=0; “Tempo de reparo em To 


neta=514.82; % Fator de escala 


beta=2.264; ZFator de forma 


gama=0.8858; “Função gama para beta igual a 2,264 
alfa =0.592; “Parâmetro lognormal 
betalog=8.631; “Parâmetro lognorma 


while rel<Ts 


1f To<cont<Tot 


r=rand (1); 


tef= log(l+r)* (neta*gama); TEF=tef+To; 


1f TEF>Tot 


TT=Tot; 


TR=r*betalog*exp (alfa”2/2);TR1=0; 


elseif To<TEF<Tot TR=r*betalog*exp (alfa”“2/2); TRI=TR; 


TTI=TEF+TRI; 


end cont=[cont, TT]; 


Tfalha=[Tfalha, TEF]; 


E 


[reparo=[Treparo, TR] 


; Tpl=0.8*Treparo; 


end Tp=0.8*TR; 


To=To+(228+Tp); 


Tot=Tot+(228+Tp); 


rel=rel+(228+Tp); 


end 


Com isso, obteve-se valores para Tfalha (O tempo em que ocorreram as falhas) e 
Treparo (Tempo de duração do reparo) em horas, expressos em ordem crescente (De cima 


para baixo, da esquerda para a direita) nas tabelas 21 e 22, respectivamente. 


ER RR E 2 


1956,86 1209760] 22583,00]  32801,22 
1971,08[ 1238580] 2275200] 352163 
202425] 12882,16] 2552048) 03050,8 
219349] 1302547] 2552280] 35644,47 
2646,74]  19108,75] 25560,05]  34161,92 
2900,40/  13150,81]  23695,18] G4211,0 
298376] 13664,98] 2405581] 94551,17 
345020 15755,05| 24155,97] 94716,54 
3499,14] 14117,00] 2455441] 94098,17 
4099,46] 1420758] 2461092] 95010,0 
4156,76]  14645,08] 2498482] 55171,92 
4176,32]  14715,98] 24988,00] 35360,14 
442606] 1483830] 2519521] 055912 
4512,55] 1494016 25390,98| — 35800,7 
497126] 1509961] 2576942] 36121,45 
5309,88] 1561827] 2618950] 36327,25 
5467,16] 1613201] 2632209] 36441,17 
572597] 1634051] 26995,18] 668826 
5880,95] 16361,16] 2672032] 37105,15 
606484] 1648251] 269292] 572268 
6304,47] 1684006] 2697767] 375593 
673402] 1694642] 27120,07]  37888,47 
651542] 1720681] 2786925] 680095 
724550] 1741490] 2704840 38052, 
735555] 1780276] 27980,08]  38495,41 
705604] 18144,43]  28963,45] 38567. 
B042,72]  18329,19] 2860056] 38678,0 
B918,44] 1833830] 288037] 350658 
6456,68 1864753] 2899749] 99712,97 
856822] 1689837] 2917265] 3973472 
629,40] 18940,35] 2973672] 95898,0 
809,11] 1927581] 29780,26| 398814 
907708] 1951963] 29965,45]  40785,04 
969,76] 1952548] — 00058,15 
9672,82]  20099,00[ 3034157 
9957,26] 20145,78] 3040768 | 
[56774] 2085186] arod4 26) | 
[1017967] 2077509] 312067) 


Tabela 27 - Tempo de falha (em horas) gerado por Monte Carlos, em ordem 
crescente. 
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Tabela 28 - Tempo de reparo (em horas) gerado por Monte Carlos, em ordem 
crescente. 


9 - ANALISE DA AMOSTRA DOS DADPS CENSURADOS 


A presença de censura configura um conjunto incompleto de dados. Isso 
porquealguns dados de falha ou reparo podem ser perdidos durante a amostragem, 
ou mesmo ocorrer uma censura à direita, resultado de paradas em testes de falha. 
No caso da máquina hipotética deste relatório, foi assumido que os dados 
censurados surgiram devido à primeira causa citada. Nesse sentido, a maioria das 
análises feitas anteriormente para dados completos foi desenvolvida novamente 
para seis dados censurados, mostrados na tabela 23, que abrange tempos de falha 


e de reparo. Os dados censurados foram marcadoscom tarja vermelha. 


1 20:03 431.59 04:39 
2 00:44 455.64 04:03 
3 20:55 503.28 30:19:00 
4 06:19 504.19 09:32 
5 09:13 504.21 09:57 
6 19:05 528.25 05:12 
7 14:19 528.53 15:02 
8 16:44 552.53 13:24 
9 03:48 576.18 08:43 
10 06:42 597.36 06:22 
11 21:34 624.61 21:38 
12 14:24 671.61 06:22 
13 11:34 671.79 09:43 
14 12:43 671.94 06:58 
15 12:30 672.29 01:40 
16 05:17 695.71 10:01 
17 07:51 695.86 10:27 
18 00:54 718.32 15:54 
19 21:23 718.98 08:53 
20 21:18 719.64 02:36 


21 08:12 719.88 08:56 
22 16:14 72041 04:25 
23 03:31 746.57 14:20 
24 04:17 768.18 11:58 
25 08:28 768.49 0.25 
26 11:59 768.84 09:01 
27 17:59 791.30 08:57 
28 10:50 792.07 11:14 
29 07:17 814.99 05:58 
30 10:02 816.10 07:04 
31 18:05 816.31 11:33 
32 08:26 839.91 00:27 
33 01:19 839.91 10:46 
34 03:51 840.18 06:26 
35 17:32 840.49 18:42 


Tabela 28 - Representação dos dados censurados (Tempos em horas). 


9.1 Método de Ajuste Rank 


Para dados censurados, os métodos de rank médio e mediano não são 
eficazes para a estimativa das características das funções de falha, uma vez que 
a presença dessesdados, além de não fornecerem a informação desejada, põe em 
conflito o posicionamentodos demais dados, quando ordenados. Nesse sentido, 
o método do ajuste de rank torna-seo método empírico mais apropriado para 
estes casos. O método consiste em ajustar os índices dos dados ordenados, 


utilizando as seguintes equações: 


(n + 1) — lt;- 


Incrementoderamk="—————————>—>——>—— > 
1 + númerodeunidadesalémdacensurada 


it, = it,-1 + Incrementoderank 


A tabela 28 foi desenvolvida a partir da lógica anterior. 


431.59 


| 04:39 0 | 
4 504.19 09:32 1.06061 2.06061 
5 504.21 09:57 0 3.06061 
6 528.25 05:12 0 4.06061 
7 528.53 15:02 0 5.06061 
8 552.53 13:24 0 6.06061 
9 576.18 08:43 0 7.06061 
10 597.36 06:22 0 8.06061 
1 624.61 21:38 0 9.06061 
12 
13 671.79 09:43 1.12247 10.18308 
14 671.94 06:58 0 11.18308 
15 672.29 01:40 0 12.18308 
16 695.71 10:01 0 13.18308 
17 695.86 10:27 0 14.18308 
18 
19 718.98 08:53 1.21205 1599518 
20 719.64 02:36 0 16.39513 
21 
200 720.41 04:25 1.30699 17.70212 
23 746.57 14:20 0 18.70212 
24 
25 768.49 0.25 1.44149 20.14361 
26 768.84 09:01 0 21.14361 
2] 791.30 08:57 0 22.14361 
28 792.07 11:14 0 23.14361 
29 814.99 05:58 0 24.14361 
30 816.10 07:04 0 25.14361 
31 816.31 11:33 0 26.14361 
2, 839.91 00:27 0 27.14361 
33 839.91 10:46 0 28.14361 


34 


840.18 


06:26 


29.14361 


ao 


840.49 


18:42 


30.14361 


Utilizando os novos índices, as funções relacionadas a falha e reparo foram 


novamente calculadas e plotadas (tabelas 28 a 30, ). 
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Analise de reparo 


Taxa de 
ti Manutebalidade acumulativa | Densidade Falha 
1 2.1 0.006448291 0.016425702 | 0.016532307 
4|3.5333333 0.067252738 0.070126752 | 0.075183016 
5 | 4.1666667 0.11843574 0.090534448 | 0.102697502 
6 | 4.8333333 0.183946592 0.104740379 | 0.128349907 
7|5.2833333 0.232378444 0.110017659 | 0.143322784 
8 5.65 0.273134393 0.112025627 | 0.154121512 
9 5.75 0.284346627 0.112254483 | 0.15685594 
10 6.4 0.357112022 0.110947826 | 0.17257723 
11|6.5333333 0.371851227 0.110171518 | 0.175390804 
13| 6.9666667 0.418881297 0.106730378 | 0.183663642 
14 205) 0.448717906 0.103865419 | 0.188407024 
15 | 7.3166667 0.45561627 0.103134715 | 0.189452237 
16 7.6 0.484369625 0.099832031 | 0.193611618 
17 | 8.2833333 0.549594771 0.090957514 | 0.201945954 
19 8.6 0.577704724 0.086611286 | 0.205096506 
20 | 8.7833333 0.593346769 0.084071099 | 0.206739041 
22 8.95 0.607162433 0.081759213 | 0.20812473 
25 7/5) 0.668160044 0.070852188 | 0.213513132 
26 9.85 0.675177281 0.069530122 | 0.214055601 
27 | 10.583333 0.722704309 0.06025267 | 0.217286717 
28| 11.583333 0.777176855 0.049008188 | 0.219942087 
29 11.65 0.780420286 0.048320707 | 0.22005998 
30 | 12.333333 0.811131304 0.041723868| 0.22091468 
31 | 13.166667 0.84290866 0.034751112 | 0.221215964 
32 13.95 0.86788525 0.029186179 | 0.22091537 
33 | 17.316667 0.936721358 0.013648357 | 0.215686628 
34 18.55 0.951415322 0.010342233 | 0.21287026 
35 | 19.983333 0.964097981 0.0075147 | 0.20931136 


Tabela 29 - Funções dos dados de reparo (censurados) 


Taxa de Falha(reparo) 


Manutebalidade acumulativa 


0 5 10 15 


— Taxa de Falha 


Manutebalidade acumulativa Densidade 


FIGURA 36 a 39: Gráfico Densidade de Reparo 


Observa-se que é possível manter a integridade das funções, sem o 
aparecimentode “buracos” que surgiriam no eixo das abcissas devido à censura. 
O comportamento das curvas é bem semelhante ao obtido nos ranks médio e 
mediano, uma vez que ainda trata-se de uma estimativa inicial (empírica). 

O conjunto agora censurado totaliza 29 dados. As mesmas equações 
utilizadas para o cálculo de TMEF e TMPR na seção? foram utilizadas, para a 


estimativa da disponibilidade. O resultado obtido encontra-se na tabela 2/77. 


Tabela 30 Disponibilidade estimada para dados censurados com 
ajuste de rank 
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